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Abstract 

New  explosive  compounds  that  are  less  sensitive  to  shock  and  high  tem¬ 
peratures  are  being  tested  as  replacements  for  TNT  and  RDX.  Two  of  these 
explosives,  DNAN  (2,4-dinitroanisole)  and  NTO  (3-nitro-i,2,4-triazol-5- 
one),  have  good  detonation  characteristics  and  are  the  main  ingredients  in 
a  suite  of  insensitive  munitions  (IM)  explosives  that  are  being,  or  soon  will 
be,  fielded.  Both  compounds,  however,  are  more  soluble  than  either  TNT 
or  RDX;  and  research  has  shown  that  both  have  some  human  and  envi¬ 
ronmental  toxicity.  Data  on  their  fate  is  needed  to  determine  if  DNAN  and 
NTO  have  the  potential  to  reach  groundwater  and  be  transported  off  base, 
an  outcome  that  could  create  future  contamination  problems  on  military 
training  ranges  and  trigger  regulatory  action.  In  this  study,  we  measured 
how  quickly  DNAN  and  NTO  dissolve  in  IM  formulations  and  how  solu¬ 
tions  of  these  IM  explosives  interact  with  different  types  of  soils.  Because 
both  dissolution  and  solution-soil  interactions  are  determined  by  a  suite 
of  parameters,  we  are  using  a  multifaceted  approach  to  studying  these 
processes.  Given  a  mass  of  IM  compounds  scattered  on  the  range,  our 
work  will  help  determine  the  dissolved  IM  masses,  their  subsequent 
transport  and  fate,  and  their  likelihood  of  reaching  groundwater. 
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1  Introduction 

1.1  Background 

Transforming  solid  high  explosives,  scattered  onto  range  soils  as  a  result  of 
incomplete  detonations,  into  solutions  of  explosive  components  that  reach 
groundwater  involves  a  series  of  steps.  These  steps  include  dissolution, 
phototransformation  (both  of  the  solid  and  of  the  aqueous  solutions 
formed  by  precipitation  interacting  with  the  solid),  and  complex  interac¬ 
tions  of  the  dissolved  explosives  with  soil  constituents  during  transport 
from  the  surface  to  the  groundwater.  If  these  compounds  and  their  degra¬ 
dation  products  are  innocuous,  their  transport  to  groundwater  might  not 
matter.  Unfortunately,  toxicology  data  for  DNAN  (2,4-dinitroanisole)  and 
NTO  (3-nitro-i,2,4-triazol-5-one)  show  that  while  NTO  has  low  mammali¬ 
an  toxicity  (London  and  Smith  1985),  DNAN  is  more  toxic  to  mammals 
than  TNT  (2,4,6-trinitrotoluene)  (Davies  and  Provatas  2006),  can  inhibit 
seed  germination  and  plant  growth,  and  even  is  used  as  a  pesticide 
(Dumitras-Hutanu  et  al.  2009).  Dodard  et  al.  (2013)  also  showed  that  it  is 
toxic  to  bacteria  and  earthworms.  Furthermore,  both  DNAN  and  NTO  can 
form  toxic  transformation  products  (Le  Campion  et  al.  1999b;  Davies  and 
Provatas  2006). 

DNAN  and  NTO  are  readily  soluble  in  water,  and  their  solubility  is  consid¬ 
erably  higher  than  those  of  TNT  and  RDX  (i,3,5-hexahydro-i,3,5-trinitro- 
1,3,5-triazine)  (Table  1;  Figure  1).  When  these  compounds  are  part  of  a 
formulation,  their  dissolution  will  depend  not  only  on  their  individual  sol¬ 
ubility  and  dissolution  rates  but  also  on  the  fraction  of  each  component 
exposed  to  water  (Lever  et  al.  2005;  Dontsova  et  al.  2006;  Taylor  et  al. 
2009a,  2009b).  In  Composition  B  (Comp  B— 60%  RDX,  39%  TNT,  1% 
wax),  for  example,  the  more  soluble  TNT  dissolves  from  the  surface,  ex¬ 
posing  progressively  more  of  the  slowly  dissolving  RDX  crystals  (Lever  et 
al.  2005).  Like  Comp  B,  the  insensitive  munitions  (IM)  formulations  are 
produced  by  melt  casting  where  one  or  more  explosive  fillers  (such  as  AP 
[ammonium  perchlorate],  RDX,  NTO,  or  NQ  [nitroguanidine])  are  mixed 
with  a  melted  energetic  binder  (TNT  or  DNAN)  (Pelletier  et  al.  2010; 
Rutkowski  et  al.  2010).  However,  in  IMs,  the  filler,  NTO,  has  a  higher  sol¬ 
ubility  than  the  binder,  DNAN.  Depending  on  the  amount  of  each  com¬ 
pound  in  the  formulation,  an  initial  loss  of  NTO  could  be  followed  by  dis¬ 
solution  of  both  components;  or  DNAN  may  protect  NTO  from  dissolving. 
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Solubility,  along  with  particle  size  and  climatic  conditions  (rainfall  and 
temperature),  controls  the  dissolution  rate  of  solid  explosives  (Lever  et  al. 
2005;  Taylor  et  al.  2009a,  2009b).  Our  aim  is  to  parameterize  what  con¬ 
trols  dissolution  of  IM  formulations. 


Table  1.  Selected  environmentally  relevant  chemical  and  physical  properties  of  DNAN,  NTO, 
TNT,  and  RDX,  including  solubility  at  25°C,  octanol-water  partition  coefficient  (, Kow ),  and  soil 
organic  carbon  adsorption  coefficient  ( Koc ). 


Property 

DNAN 

NTO 

TNT 

RDX 

Solubility,  mg  L_1 

276. 2a  1 

16642. 0C  2 

100. 5a  3 

59.9a3 

Log  Kow 

1.7-1.92b,  1.64a4 

0.37-1.03b4 

1.6-1.84  5 

0.81-0.87  5 

Log  Koc 

3.11a 6,  2.2b7 

3.03a  6,  2.1b7 

3.2 3 

0.88-2.4 3 

a  measured; b  estimated; c  interpolated  from  measured  values 

1(Boddu  etal.  2008); 2 (Spear  etal.  1989); 3 (Brannon  and  Pennington  2002);  4(Sokkalingam  etal.  2008); 5 (Yoon  etal.. 
2005); 6  (R.  Pesce-Rodriguez,  unpublished); 7 (Chakka  etal.  2010). 


Little  is  known  about  the  phototransformation  of  DNAN  and  NTO  either 
as  solids  or  in  solution.  This  process  is  important  because  sunlight  might 
chemically  alter  the  surface  of  the  IM  compound  to  produce  other  com¬ 
pounds  with  different  solubilities  or  toxicities.  Taylor  et  al.  (2010)  discov¬ 
ered  this  to  be  the  case  for  TNT  and  RDX.  Preliminary 
phototransformation  studies  of  DNAN  and  NTO  solutions  made  using  an 
accelerated  program  in  an  Atlas  05000  system  (uses  a  xenon  lamp  and  a 
6-month  cycle  equivalent  to  2  years  of  real-time  exposure  in  southern 
Florida)  indicated  a  half-life  of  less  than  3  days  for  DNAN  and  NTO  (R. 
Pesce-Rodriguez,  unpublished  data).  These  data  suggest  that  in  the  pres¬ 
ence  of  light,  DNAN  and  NTO  in  solution  are  more  stable  than  TNT  but  are 
less  stable  than  RDX.  Rao  et  al.  (2013)  observed  even  faster 
phototransformation  for  DNAN,  with  half-lives  ranging  between  0.11  and 
1.51  days.  Hydrolysis  studies  found  that  without  sunlight,  DNAN  and  NTO 
are  stable  in  neutral,  acidic,  and  basic  solutions  (R.  Pesce-Rodriguez,  un¬ 
published  data).  Observed  products  of  phototransformation,  according  to 
Le  Champion  et  al.  (1999a),  include  nitrites,  nitrates,  and  carbon  dioxide. 
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Figure  1.  Plots  of  the  solubility  of  (a)  DNAN  and  (b)  NTO  as  a  function  of  temperature.  Data 
are  from  studies  conducted  at  the  U.S.  Army  Construction  Engineering  Research  Laboratory 
(CERL)  (Boddu  et  al.  2008),  the  Army  Research  Laboratory  (ARL)  (R.  Pesce-Rodriguez, 
unpublished),  the  Defense  Science  and  Technology  Organization  (Smith  and  Cliff  1999),  and 

others. 
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Outdoor  dissolution  tests  of  Comp  B  and  TNT  showed  that 
phototransformation  of  the  explosive  pieces  was  a  major  loss  pathway 
(Taylor  et  al.  2010).  One-half  to  two-thirds  of  the  mass  dissolved  was  un¬ 
accounted  for,  and  losses  scaled  closely  with  surface  area  (Figure  2).  A 


ERDC/CRREL  TR-14-23 


4 


similar  phototransformation  process  may  be  occurring  on  the  IM  pieces 
we  placed  outside  to  dissolve. 

For  TNT  and  Comp  B,  the  dissolved  colored  surface  products  were  associ¬ 
ated  with  large  pre-solvent  peaks.  Analyses  of  the  water  samples  showed 
the  presence  of  2-amino-4,6-dinitrobenzoic  acid  at  concentrations  similar 
to  or  larger  than  the  TNT  concentrations.  As  2-amino-4,6-dinitrobenzoic 
acid  is  a  known  transformation  product  of  TNT,  this  compound  could  ac¬ 
count  for  the  missing  TNT  mass.  Because  2-amino-4,6-dinitrobenzoic  acid 
is  very  polar,  it  is  not  quantified  by  Method  8300B  (U.S.  EPA  2006)  and 
elutes  in  the  pre-solvent  peak  area  of  the  chromatogram  (Taylor  et  al. 
2010). 

We  proposed  identifying  and  quantifying  any  phototransformation  prod¬ 
ucts  occurring  on  IM  chunks,  but  this  portion  of  our  project  was  not  fund¬ 
ed.  If  monies  are  available,  we  think  that  it  is  important  to  do  this  work  as 
these  phototransformation  products  are  soluble  and  mobile;  and  we  do  not 
know  what  they  are. 


Figure  2.  Mass  loss  measured  by  high  precision  liquid  chromatography  (HPLC)  versus 
measured  mass  loss  by  electronic  balance  for  TNT  and  Tritonal  (a  TNT-aluminum  explosive) 

samples  (Taylor  et  al.  2010). 
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Once  in  solution,  DNAN,  NTO,  and  their  transformation  products  will  ex¬ 
perience  reactive  transport  through  the  soil.  En  route  to  ground  water, 
they  can  undergo  irreversible  and  reversible  adsorption  by  different  min¬ 
eral  and  organic  phases  in  the  soil,  transformation,  volatilization,  and  bio¬ 
uptake.  The  importance  of  each  of  these  processes  for  the  fate  of  a  com¬ 
pound  can  be  evaluated  using  the  soil  adsorption  coefficient  (Kd),  Koc,  Kow, 
Henry’s  Law  constant,  the  one-electron  standard  reduction  potential,  and 
the  first-order  transformation  rate  constant  (k),  among  others.  Table  l 
provides  values  for  some  of  these  parameters  for  DNAN  and  NTO  and  for 
common  explosives  TNT  and  RDX. 

The  affinity  of  the  energetics  to  be  sorbed  onto  soil  surfaces  is  character¬ 
ized  by  the  soil  adsorption,  or  partitioning,  coefficient,  Kd.  In  particular, 
organic  matter  (OM),  phyllosilicate  clays,  and  iron  and  aluminum  oxides 
and  hydroxides  adsorb  organic  compounds  due  to  their  high  surface  areas 
(Dontsova  et  al.  2009a).  As  values  of  Kd  increase,  the  chemical  resides 
mostly  on  soil  surfaces  and  little  is  transported  downward  in  the  moving 
pore  water.  As  most  of  the  adsorption  of  organic  contaminants  in  soils  can 
be  attributed  to  soil  organic  carbon,  Kd  values  are  often  normalized  to  soil 
organic  carbon  content,  and  the  resulting  parameter  is  called  the  soil  or¬ 
ganic  carbon  adsorption  coefficient  (Koc).  This  parameter  is  estimated  for 
the  compound  and  then  used  to  calculate  adsorption  to  other  soils,  de¬ 
pending  on  their  carbon  content.  There  is  little  information  about  soil  ad¬ 
sorption  coefficients  for  the  IM  compounds.  Recent  work  reports  Kd  of 
2.7  cm3  g-i  for  DNAN  and  less  than  0.1  cm3  g-i  for  NTO  in  soil  (2.08%  or¬ 
ganic  carbon  [OC]),  indicating  that  natural  attenuation  maybe  an  im¬ 
portant  mechanism  for  DNAN  remediation  (Hawari  et  al.  2011,  2012). 
However,  more  work  is  needed  to  measure  how  chemical  and  biological 
transformation  would  affect  the  measured  Koc  values  (as  in  the  standard 
ASTM  E1195-01  method  [ASTM  International  2008],  adsorption  is  deter¬ 
mined  by  the  difference  between  initial  and  equilibrium  concentrations  of 
a  compound  in  solution)  and  to  determine  the  mechanisms  involved.  For 
example,  RDX  tends  to  partition  to  OC  (Tucker  et  al.  2002;  Dontsova  et  al. 
2009a)  while  TNT  exhibits  very  complex  sorption  behaviors,  which  in¬ 
clude  adsorption  in  clay  interlayer  surfaces  (Haderlein  et  al.  1996)  and  ir¬ 
reversible  adsorption  to  OM  following  transformation  (Thorn  and 
Kennedy  2002). 

Unlike  most  other  explosives,  which  are  polar  but  non-ionic  compounds, 
NTO  is  an  acid  with  an  acid  dissociation  constant  (p Ka)  of  3.7-3.76 
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(Chipen  et  al.  1966;  Le  Campion  et  al.  1997)  and  will  be  negatively  charged 
at  environmentally  relevant  pHs  (Smith  and  Cliff  1999).  As  both  organic 
and  mineral  soil  surfaces  tend  to  have  a  net  negative  charge,  we  would 
predict  low  affinity  of  NTO  for  soils  and  greater  mobility  in  the  environ¬ 
ment.  However,  negatively  charged  organics  can  be  adsorbed  in  soils 
through  cation  bridging. 

The  octanol-water  partition  coefficient  is  a  measure  of  the  tendency  of  the 
chemical  to  bio-accumulate  (Boddu  et  al.  2008)  and  an  indicator  of  its  af¬ 
finity  for  soil  OM.  Relatively  low  reported  Kow  values  (Table  1)  indicate  that 
these  compounds  would  not  be  strongly  adsorbed  by  OM  in  soils  through 
non-polar  interactions.  However,  they  can  still  experience  specific  sorp¬ 
tion,  as  has  been  shown  for  TNT  and  the  DNTs  (2,4-  and  2,6- 
dinitrotoluene)  (Thorn  et  al.  2001;  Thorn  et  al.  2008).  Based  on  the  Kow 
values,  one  would  expect  NTO  to  be  more  mobile  than  DNAN  and  to  ex¬ 
hibit  similar  environmental  behaviors  to  RDX. 

Relatively  low  volatility,  as  measured  by  Henry’s  constants  (logH  of  -3.25 
to  -4.40  for  DNAN  and  -11.38  for  NTO)  (Sokkalingam  et  al.  2008),  indi¬ 
cates  that  transport  in  a  gas  phase  will  not  be  important  for  the  total  mass 
flux  of  IM  compounds  in  soils.  One-electron  standard  reduction  potential 
(Em)  is  a  measure  of  how  easily  nitro  groups  in  these  compounds  reduce  to 
amino  groups  (Uchimiya  et  al.  2010),  a  common  first  step  in  transfor¬ 
mation  of  many  explosive  compounds.  Boparai  et  al.  (2010)  observed  fast¬ 
er  reduction  for  explosives  having  more  positive  one-electron  reduction 
potentials.  Measured  electron  standard  reduction  potential  for  DNAN 
(-0.40  volt)  was  similar  to  2,4-DNT,  which  is  readily  reduced  in  soil  envi¬ 
ronments.  Note  that  transformation  products  of  the  compound  can  have 
different  environmental  and  toxicity  properties  than  the  original  contami¬ 
nant. 

Microbial  degradation  studies  of  NTO  in  soils  show  nitroreduction  fol¬ 
lowed  by  the  ring  cleavage  of  the  primary  amine,  5-amino-i,2,4-triazol-3- 
one  (ATO)  (Le  Campion  et  al.  1999b).  At  optimal  conditions,  the 
nitroreduction  of  NTO  was  complete  in  24  hr,  while  the  degradation  of 
ATO  required  2  weeks  of  incubation.  The  end  products  of  the  biodegrada¬ 
tion  were  carbon  dioxide  (40%),  urea,  and  a  polar  compound  assumed  to 
be  hydroxyurea. 
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In  summary,  there  is  insufficient  information  about  soil  adsorption  for  the 
IM  compounds.  Both  the  extent  of  sorption  and  the  mechanisms  involved 
are  not  yet  known.  Kow  has  been  estimated  and  measured  (Sokkalingam  et 
al.  2008),  but  the  Koc  values  that  are  currently  available  were  computa¬ 
tionally  derived  (Chakka  et  al.  2010)  from  Kow,  which  can  result  in  unreli¬ 
able  values,  particularly  for  polar  compounds;  and  soil  adsorption  was 
measured  for  only  one  soil  (Hawari  et  al.  2011,  2012),  making  it  difficult  to 
generalize  how  IM  compounds  will  behave  in  other  soils.  To  address  these 
knowledge  gaps,  we  will  measure  DNAN  and  NTO  adsorption  coefficients 
( Kd )  for  multiple  soils  having  a  range  of  properties  and  will  normalize  the 
values  to  their  organic  carbon  content  to  determine  Koc.  If  these  parame¬ 
ters  are  known,  the  fate  of  the  chemical  in  the  soil  system  can  be  modeled. 
Because  soils  are  an  important  mechanism  for  natural  remediation  of  or¬ 
ganic  and  inorganic  contaminants,  values  for  Kd,  k,  and  Koc  are  critical  for 
understanding  the  fate  and  transport  of  the  IM  compounds. 

1.2  Objectives  and  approach 

We  seek  to  measure  how  quickly  DNAN  and  NTO  dissolve  in  IM  formula¬ 
tions  and  how  solutions  of  these  IM  explosives  interact  with  different  types 
of  soils.  Because  both  dissolution  and  solution-soil  interactions  are  de¬ 
termined  by  a  suite  of  parameters,  we  are  using  a  multifaceted  approach  to 
study  of  these  processes.  Given  a  mass  of  IM  compounds  scattered  on  the 
range,  our  work  will  help  determine  the  dissolved  IM  masses,  their  subse¬ 
quent  transport  and  fate,  and  their  likelihood  of  reaching  groundwater. 

To  measure  the  dissolution,  phototransformation,  and  soil  adsorption 
properties  of  DNAN,  NTO,  and  IM  formulations  that  contain  them 
(IMX101,  IMX104,  and  PAX21),  we 

1.  used  light  and  electron  microscopy  to  characterize  particles  of  DNAN, 
NTO,  and  detonation  residues  of  IM  formulations  that  contain  them; 

2.  measured  the  dissolution  of  detonation  residues  of  IM  formulations  in  a 
controlled  laboratory  setting  and  in  the  field  where  they  were  exposed  to 
rainfall  and  sunlight; 

3.  measured  the  phototransformation  of  DNAN  and  NTO  in  solution  and 
from  solid  IM  formulations  set  outside  to  degrade  and  dissolve; 

4.  determined  adsorption  behavior  of  IM  compounds  in  a  range  of  soils;  and 

5.  measured  the  transport  behavior  of  DNAN  and  NTO  in  representative 
soils,  which  helped  us  correlate  their  transport  and  fate  behavior  to  soil 
properties. 
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2  Dissolution  Studies 

2.1  Materials  and  methods 

2.1.1  IM  particle  collection 

All  particles  came  from  partially  detonated  rounds  and  were  collected  in 
the  field.  We  obtained  detonation  residues  of  IMXioi  from  tests  conduct¬ 
ed  on  105  mm  howitzer  rounds  by  National  Technical  Systems  and  then 
collected  PAX21  and  IMX104  residues  after  low-order  detonations  of  60 
and  81  mm  mortars,  respectively  (Walsh  et  al.  2013,  2014). 


Table  2.  IM  formulations  containing  DNAN  and  NTO  in  advanced  stages  of  development  or 

ready  for  fielding. 


IM 

Replaces 

Components 

Used  In 

IMXIOI 

TNT 

43.5%  DNAN,  36.8%  NQ, 
19.7%  NTO 

Qualified  by  the  U.S.  Army  as  the  main  fill 
in  the  155  mm  M795  Artillery  Projectile 

IMX104 

Comp  B 

53%  NTO,  31.7%  DNAN, 
15.3%  RDX 

Used  by  the  U.S.  Army  for  60,  81,  and  120 
mm  mortars  and  sub-munitions 

PAX21 

Comp  B 

34%  DNAN,  30%  AP, 
30.75%  RDX 

Qualified  and  fielded  in  the  U.S.  Army  60 
mm  mortar  system 

Some  data  from  Fung  et  al.  (2010),  Wilson  (2007),  Pelletier  et  al.  (2010),  and  Coppola  (2007). 


2.1.2  Imaging 

We  used  a  SkyScan  1172  X-ray  microtomograph  to  image  the  IM  formula¬ 
tions.  We  ran  the  instrument  at  a  40  kV  by  using  a  250  pA  source  current. 
A  1.3  megapixel  X-ray  camera  took  a  sequence  of  shadow  transmission  im¬ 
ages  (radiographs)  of  each  IM  particle  as  it  rotated  on  the  stage.  From  the¬ 
se  radiographs,  cross  sectional  images  perpendicular  to  the  z-axis  of  the 
particle  were  reconstructed  using  a  modified  Feldkamp  cone-beam  algo¬ 
rithm,  creating  a  complete  three-dimensional  (3-D)  representation  of  the 
internal  microstructure  of  the  particle.  We  quantified  the  areas  of  the  ma¬ 
trix  and  of  the  crystal  components  in  each  section  using  Image  J  (Rasband 
2012).  We  calculated  the  volumes  by  multiplying  these  area  measurements 
by  the  spacing  between  cross  sections  (4  to  7  pm).  This  technique  required 
no  sample  preparation. 

2.1.3  Indoor  drip  tests 

We  used  well-controlled  laboratory  drip  experiments  to  obtain  dissolution 
data  on  particles  of  detonated  IMX101,  IMX104,  and  PAX21.  Our  experi- 
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mental  design  mimicked  field  conditions  on  training  ranges,  where  rainfall 
dissolves  spatially  isolated  particles  and  chunks  of  explosives  scattered  on 
the  soil  surface  (Taylor  et  al.  2009a).  We  placed  a  massed  piece  of  the  ex¬ 
plosive  on  a  glass  frit,  dripped  water  on  the  particle  at  a  given  rate,  and 
analyzed  the  effluent  using  HPLC  (Figure  3a).  No  soil  was  involved  in  the 
tests,  allowing  us  to  determine  the  dissolution  as  a  function  of  the  particle 
mass  and  the  water  volume  interacting  with  the  particle.  Figure  3a  shows 
the  syringe  pump  used  in  these  tests. 


Figure  3.  Photos  of  (a)  the  laboratory  drip  tests  and  of  (b)  the  outdoor  setup. 


2.1.4  Outdoor  drip  tests 

We  also  conducted  dissolution  tests  using  centimeter-size  chunks  of  high 
explosives  that  we  placed  outside  to  weather  under  natural  conditions  at 
the  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL)  facility 
in  Hanover,  NH  (Figure  3b).  These  tests  closely  simulated  the  dissolution 
of  isolated  high  explosive  pieces  on  range  soils. 

In  June  2013,  we  placed  5  chunks  of  IMX101,  5  chunks  of  IMX104  and  2 
chunks  of  PAX21  outside  in  4  cm  diameter  Buchner  funnels  (Figure  3b). 
The  funnels  were  attached  to  1 L  bottles  with  a  #4  rubber  stopper  that  was 
fitted  with  two  holes— one  for  the  funnel  stem  and  the  other  for  air  ex¬ 
change.  The  bottles  fitted  snuggly  into  an  insulated  wooden  box.  The  box 
kept  the  bottles  from  tipping  over  or  being  damaged  and  also  helped  to 
keep  them  cool  in  summer  and  warm  in  winter. 

Rainwater  or  snowmelt  that  interacted  with  the  high  explosives  was  col¬ 
lected  in  the  bottles.  Periodically,  we  exchanged  the  bottles  for  clean  ones 
and  measured  the  volume  and  high-explosive  concentration  in  the  water. 
During  the  winter  months,  we  placed  the  bottles  in  plastic  bags  to  save  the 
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sample  even  if  the  bottle  should  break  from  ice  formation.  During  each 
collection,  we  photographed  the  pieces  of  high  explosives  in  situ  to  docu¬ 
ment  changes  in  their  appearance  and  size.  This  setup  allowed  the  explo¬ 
sives  to  be  exposed  to  conditions  similar  to  those  experienced  by  explo¬ 
sives  on  a  range— where  rain,  snow,  sun,  and  freeze-thaw  cycles  weather 
them— while  allowing  us  to  collect  and  analyze  the  dissolved  high  explo¬ 
sives  and  to  monitor  changes  in  their  appearance.  We  set  up  our  experi¬ 
ment  in  an  enclosed,  locked  site  next  to  a  dedicated  rain  gauge  and  near 
the  CRREL  weather  station,  which  records  rainfall  totals,  temperature, 
and  solar  irradiance.  We  downloaded  the  rain  data  monthly. 

2.1.5  Analytical  methods 

We  analyzed  by  HPLC  all  solutions  from  the  drip  and  the  outdoor  tests 
and,  using  an  ion  selective  probe,  analyzed  solutions  from  PAX21  particles 
for  perchlorate.  We  used  two  HPLC  methods:  Method  8330  (U.S.  EPA 
2006)  to  quantify  DNAN,  RDX,  and  HMX  (a  manufacturing  impurity  of 
RDX)  and  a  second,  described  below,  to  quantify  NTO  and  NQ.  For  both 
methods,  we  ran  a  1  or  10  ppm  standard  every  ten  samples  and  ran  blanks 
before  and  after  each  standard  run  to  prevent  carryover. 

We  used  the  U.S.  Environmental  Protection  Agency’s  SW-846  Method 
8330  (U.S.  EPA  2006)  to  quantify  DNAN,  RDX,  HMX,  their  co¬ 
contaminants,  and  breakdown  products.  We  added  DNAN  to  the  commer¬ 
cially  available  8330  standard  (Restek)  and  used  this  solution  for  calibra¬ 
tion.  We  used  a  Waters  NovaPak  C8  column  (3.9  x  150  mm)  eluted  at  1.4 
mL  min-1  (28°C)  with  85:15  water:isopropanol  mix  and  detected  by  UV  at 
295  nm.  We  added  1  mL  of  the  effluent  to  2  mL  of  deionized  water  and  to  1 
mL  of  acetonitrile  and  then  filtered  it  through  a  0.45  pm  Millipore  syringe 
filter. 

We  prepared  the  samples  for  NTO  and  NQ  by  adding  1  mL  of  the  effluent 
to  3  mL  of  acetonitrile  containing  0.1%  trifluoroacetic  acid  (TFA)  and  then 
filtered  it  through  a  0.45  pm  Millipore  syringe  filter.  The  standard  con¬ 
tained  1  mL  of  a  1  mg  L_1  solution  of  NTO  in  acetonitrile  (AcN),  1  mL  of  a  1 
mg  L_1  solution  of  NQ  dissolved  in  water,  and  2  mL  of  AcN  containing 
0.1%  TFA.  We  used  a  Thermo  Scientific  Hypercarb  column  (4.6  x  150  mm) 
eluted  at  1.4  mL  min-1  (28°C)  with  95:5  methanokwater  mix  and  detected 
by  UV  at  300  nm. 
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We  quantified  the  AP  using  a  Cole  Parmer  ion-selective  electrode.  To  gen¬ 
erate  a  standard  line  from  which  we  could  determine  the  AP  concentration 
of  the  samples,  we  prepared  and  analyzed  a  series  of  concentrations  from  a 
perchlorate  standard.  Using  a  Mettler  Toledo  SevenEasy  meter  and  a  pH 
probe,  we  measured  the  pH  of  the  water  samples;  and  we  completed  a 
three-point  calibration  by  using  pH  4,  7,  and  10  solutions. 

2.2  Results 

2.2.1  Characterized  particles  of  DNAN,  NTO,  and  detonation  residues  of  IM 
formulations 

Figure  4  shows  optical  images  of  DNAN,  NTO,  and  the  IM  formulations 
and  the  compositions  of  the  formulations.  DNAN  is  a  cream-colored  parti¬ 
cle  composed  of  many  small  crystals.  It  abrades  easily  as  evidenced  by  the 
powdered  DNAN  on  its  surface.  NTO  is  a  translucent  crystal  that  occurs  as 
about  500  pm,  equidimensional  grains,  most  of  which  look  twinned  and 
are  highly  reflective  in  visible  light.  IMX101  and  104  are  white  crystalline 
substances  with  visible  NTO  crystals  embedded  in  the  DNAN  matrix. 
PAX21  has  a  yellow  DNAN  matrix  caused  by  the  addition  of  less  than  1 
wt.%  methyl  nitroaniline  with  clearly  visible  200-400  pm  AP  crystals. 

IMX101,  IMX104,  and  PAX21  are  made  by  melt  casting,  a  process  where 
crystalline  constituents  (e.g.,  AP,  NTO,  NQ,  and  RDX)  are  added  to  a  mol¬ 
ten  matrix  (DNAN)  and  then  cooled.  Melt  casting  is  also  used  to  produce 
conventional  high  explosives,  such  as  TNT  and  Comp  B.  In  the  case  of 
Comp  B,  100  pm  crystals  of  RDX  are  added  to  molten  TNT.  For  IMX101, 
360  pm  NTO  and  300  pm  NQ  crystals  are  added  to  molten  DNAN;  for 
IMX104,  approximately  300  pm  NTO  crystals  and  4  pm  RDX  crystals  are 
added  to  DNAN;  and  for  PAX21,  200-400  pm  AP  crystals  and  RDX  crys¬ 
tals  of  two  different  sizes  are  mixed  into  DNAN  (Pelletier  et  al.  2010; 
Rutkowski  et  al.  2010). 
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Figure  4.  Optical  images  of  DNAN,  NTO,  IMX101,  IMX104,  and  PAX21  and  the  compositions 

of  the  three  IM  formulations  we  studied. 


The  IM  formulations  differ  from  Comp  B  in  several  important  ways.  First, 
most  of  their  crystal  constituents  are  much  more  soluble  than  their  DNAN 
matrices,  which  is  the  opposite  of  Comp  B  where  the  RDX  is  less  soluble 
than  the  TNT  matrix.  We  surmised  (correctly)  that  as  NTO,  NQ,  and  AP 
dissolve,  they  would  leave  holes  in  the  DNAN  matrix,  thereby  increasing 
DNAN’s  surface  area  and  its  dissolution  rate  (Figure  5).  Second,  because 
the  solubility  of  NTO,  NQ,  AP  and  DNAN  differ  from  each  other  by  orders 
of  magnitude,  we  expected  these  compounds  to  dissolve  at  different  rates 
based  on  their  solubility  and  on  the  fraction  of  their  surface  exposed  to  wa¬ 
ter.  The  3-D  relationship  of  the  constituents  and  their  relative  crystal  sizes 
in  IM  formulations  is,  therefore,  important  for  understanding  and  model¬ 
ing  the  rate  at  which  IM  particles  dissolve.  Figure  6  illustrates  an  example 
of  the  effect  differences  in  solubility  can  have  on  particle  structure.  We 
found  that  the  AP  in  PAX21  particles  quickly  dissolved,  leaving  a  DNAN- 
RDX  matrix.  As  the  DNAN  dissolved,  the  small,  approximately  10  pm  RDX 
crystals  remained. 
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Figure  5.  Images  of  (a)  the  Comp  B  surface  and  of  (b)  an  IMX104  particle.  These  images 
show  that  as  the  TNT  matrix  dissolves  from  Comp  B,  the  less  soluble  RDX  crystals  jut  from  the 
surface  and  eventually  fall  off  the  particle.  For  IMX104,  on  the  other  hand,  the  DNAN  matrix  is 
less  soluble  than  the  crystal  constituents  and  dissolution  produces  hole-riddled  particles. 


Como  B  IMX104-3 


Figure  6.  PAX21  particles:  (a)  detonation  residue;  (b)  AP  completely  dissolved;  (c)  DNAN 
almost  completely  dissolved;  and  (d)  RDX  crystals  remaining. 


To  image  the  3-D  distribution  of  the  various  components  in  IMX101, 
IMX104,  and  PAX21,  we  used  micro-computed  tomography  (pCT).  Figure 
7  shows  a  radiograph  (a,  c,  e,  and  g)  and  a  reconstructed  cross  section  (b, 
d,  f,  and  h)  for  IMX101,  IMX104,  and  PAX21.  The  same  PAX21  particle  is 
shown  twice,  before  and  after  a  half  hour  submersion  in  water.  The  round- 
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ed  AP  crystals  are  clearly  visible  in  the  PAX21  radiographs  because  AP  has 
a  higher  X-ray  attenuation  that  DNAN  and  RDX.  This  is  not  the  case  for 
IMX101  (Figure  7a)  and  IMX104  (Figure  7c)  as  the  X-ray  attenuations  for 
DNAN,  NTO,  RDX,  and  NQ  are  all  similar. 

In  contrast,  the  reconstructed  cross  sections  of  the  IM  formulation  do 
show  the  constituent  crystals.  IMX101  has  square  shaped  NQ  crystals  and 
larger,  brighter  NTO  crystals  in  a  darker  DNAN  matrix  (Figure  7b).  Cracks 
in  the  DNAN  matrix  are  also  visible.  Reconstructed  cross  sections  of 
IMX104  (Figure  7d)  also  show  the  NTO  crystals  and  fractures  in  the  parti¬ 
cle. 

Figure  7.  Shown  are  radiographs  and  reconstructed  cross  sections  of  IMX101  (a  and  b), 
IMX104  (c  and  d),  and  PAX21  (e  and  f).  Also  shown  is  the  same  PAX21  particle  after  a  half- 
hour  submersion  in  water  (g  and  h).  The  scale  bar  applies  to  only  the  cross  sections.  The 
diameter  of  the  post  used  to  mount  the  samples  for  the  radiographs  was  3  mm. 


The  AP  crystals  in  reconstructed  cross  sections  of  PAX21  are  white,  large, 
rounded  objects  in  a  grey  matrix  (Figure  yf).  The  smaller  bright  phases 
(hard  to  see)  are  100  pm  RDX  crystals.  This  particle  was  placed  in  water 
for  30  min  and  then  reimaged  (Figure  7g  and  h).  Using  the  computed  cross 
sections  taken  pre-  and  post-  submersion,  we  calculated  the  volume  of  the 
particle  and  of  the  AP  crystals  within  it.  Initially  the  PAX21  particle  con¬ 
tained  110  AP  crystals  that  occupied  5  mm3,  or  24%,  of  the  21  mm3  particle 
volume.  After  a  half  hour  in  water,  63  AP  crystals  remained;  and  these  oc¬ 
cupied  2  mm3,  or  12%,  of  the  post-submersion  particle  volume  (16  mm3).  A 
decrease  in  the  percent  AP  volume  after  submersion  confirms  that  the  AP 
dissolved  preferentially.  The  AP  crystals  within  0.5  mm  from  the  periphery 
dissolved,  and  those  crystals  0.6  mm  from  the  edge  partially  dissolved. 
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The  sequence  of  cross  sections  showed  that  AP  crystals  are  about  50  pm 
apart  and  that  only  a  few  actually  touch. 

To  determine  if  the  fractures  we  saw  in  the  pCT  cross  sections  formed  dur¬ 
ing  detonation,  we  imaged  and  compared  four  manufactured  and  four  det¬ 
onated  particles  of  IMX101,  IMX104  and  PAX21  (Figure  8).  We  saw  inter¬ 
nal  cracks  in  the  detonated  particles  and  none  in  the  pCT  cross  sections  of 
the  manufactured  pieces,  indicating  that  detonation  fractured  the  parti¬ 
cles.  The  fractures  were  generally  in  the  DNAN  and  followed  the  perime¬ 
ters  of  the  AP,  NTO,  or  NQ  crystals.  A  few  of  the  fractures  extended 
through  the  entire  particle  and  could  be  followed  in  many  successive  cross 
sections. 


Figure  8.  Images  of  undetonated,  detonated,  and  wetted  IM 
particles.  Note  the  fractures  in  the  detonated  particles  and 
the  loss  of  components  in  the  wetted  particles. 


IMX101 


IMX104 


PAX21 


Undetonated  Detonated  Wetted 


2.2.2  Laboratory  tests  to  measure  the  dissolution  of  IM  detonation 
residues 

We  conducted  two  drip  dissolution  tests,  one  that  dripped  a  total  of  about 
20  mL  of  water  on  samples  that  were  imaged  using  pCT  and  the  second 
that  dissolved  the  particles  completely.  We  weighed  all  IM  samples  at  the 
beginning  of  the  tests  and  massed  the  pCT  samples  at  the  end  using  a  Met- 
tler  A230  balance  (1.0  ±  0.1  mg). 
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The  first  tests  involved  imaging  individual  IM  particles,  dripping  water  on 
them,  reimaging  them,  and  repeating  this  procedure  numerous  times.  We 
placed  each  particle  in  a  scored  plastic  pipette  tip.  The  plastic  tip  did  not 
interfere  with  the  pCT  scans,  it  held  the  particle  in  place  and  the  scoring 
allowed  us  to  place  the  samples  in  the  scanner  in  the  same  orientation 
each  time.  Although  the  pipette  tip  stabilized  the  particle,  some  of  the  par¬ 
ticles  did  move  as  they  dissolved.  The  water  dripped  on  the  particle  flowed 
out  the  tip  bottom  to  the  sample  vial. 

For  the  second  tests,  we  placed  each  particle  on  a  glass  frit  within  a  Buch¬ 
ner  funnel,  dripped  water  on  the  particle  at  0.5  mL  hr-1  using  a  syringe 
pump,  and  analyzed  the  effluent  using  HPLC.  Four  particles  of  IMX101 
and  IMX104  were  run  continuously  until  the  particle  dissolved  (month 
time  scales).  These  tests  are  similar  to  those  run  on  high  explosives  (Lever 
et  al.  2005;  Taylor  et  al  2009b).  We  did  not  run  similar  tests  on  PAX21  be¬ 
cause  its  use  was  discontinued  when  Walsh  et  al.  (2013)  found  that  AP  was 
not  consumed  during  high  order  detonations. 

For  both  of  these  tests,  we  massed  and  analyzed  the  Milli-Q  water  interact¬ 
ing  with  the  IM  particles.  One  major  difference  between  the  two  dissolu¬ 
tion  tests  was  that  particles  placed  in  the  pipette  tips  were  always  struck  by 
the  drops  whereas  those  in  the  glass  frit  were  only  struck  some  of  the  time. 
No  soil  was  present  in  either  type  of  test. 

2.2.2. 1  Dissolution  of  IM  compounds  tracked  by  pCT. 

Figure  9  shows  the  cumulative  mass  dissolved  in  the  water  dripped  on  the 
pCT-imaged  IMX101,  IMX104,  and  PAX21  particles.  For  IMX101,  more 
NTO  and  NQ  dissolved  than  DNAN;  and  more  NTO  dissolved  from  each 
particle  than  NQ  despite  there  being  more  NQ  in  the  particle.  For  IMX104, 
the  dissolved  masses  of  the  three  components  are  different  enough  to  be 
well  separated  on  the  plot  with  NTO  >  DNAN  >  RDX.  For  PAX21,  AP  dis¬ 
solved  quickly  and  completely,  followed  by  DNAN  and  RDX.  The  com¬ 
pounds  dissolved  in  accordance  with  their  solubility. 


ERDC/CRREL  TR-14-23 


17 


Figure  9.  Cumulative  mass  dissolved  in  the  water  dripped  on  the  pCT-imaged  (a)  IMX101,  (b) 

IMX104,  and  (c)  PAX21  particles. 
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The  jllCT  reconstructed  cross  sections  show  the  physical  changes  under¬ 
gone  during  dissolution  of  each  IM  formulation.  Figure  10  shows  a  cross 
section  of  each  formulation  and  how  these  changed  as  a  function  of  water 
interaction.  For  IMXioi,  as  the  NTO  and  NQ  crystals  dissolved,  they  left 
thin  bridges  of  DNAN  matrix  that  can  break  and  dissolve.  For  IMX104,  the 
dissolution  of  NTO  is  more  noticeable  as  NTO  composes  over  50%  of  the 
formulation.  Here  again,  the  DNAN  matrix  maintains  the  outline  of  the 
former  crystals.  Note  that  both  IMXioi  and  104  have  well  defined  cracks 
and  yet  dissolution  appears  to  proceed  from  the  outside  inward  and  not 
preferentially  along  the  fractures.  For  PAX21,  the  bright  AP  crystals  are 
dissolved  after  contact  with  only  a  small  volume  of  water  (2  mL).  Although 
fractures  are  present  in  PAX21,  these  not  as  well  developed  as  in  the  other 
two  formulations;  and  it  is  not  clear  how  the  water  accessed  the  AP.  After 
the  AP  dissolved,  the  RDX  crystals  could  be  seen.  They  are  the  small  grains 
that  are  slightly  brighter  than  the  DNAN  matrix.  Movies  made  from  the 
pCT  sections  are  available  in  the  supplementary  materials  of  Taylor  et  al. 
(2013,  2015). 
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Figure  10.  Change  in  the  cross  section  of  each  formulation  with  increased  water  interaction. 


From  the  |nCT  cross-section  images,  we  calculated  the  area  of  the  particle 
and  of  the  crystals  in  each  section.  We  then  estimated  the  volume  of  the 
whole  particle  and  of  its  crystal  components  by  multiplying  each  cross- 
sectional  area  by  the  spacing  between  sections  (Table  3). 

Figure  11a  plots  changes  in  particle  and  crystal  volumes  versus  water  vol¬ 
ume  for  the  four  IMX101  particles.  Note  the  large  change  in  crystal  volume 
caused  by  the  first  mL  of  water.  Over  the  course  of  the  test,  the  ratio  of  ma¬ 
trix  to  crystal  changed  from  2  to  3,  in  keeping  with  the  fact  that  NTO  and 
NQ  are  more  soluble  than  DNAN. 

We  estimated  the  density  of  the  formulations  by  plotting  the  initial  volume 
versus  mass  of  each  particle  and  fitting  a  line  to  the  data  (Figure  lib).  This 
procedure  yielded  initial  densities  of  1.69  g  cnr3  for  IMX101, 1.67  g  cnr3 
for  IMX104,  and  1.56  g  cm  3  for  PAX21.  Because  the  density  of  the  formu¬ 
lation  should  change  during  dissolution  if  dense  or  less  dense  phases  pref¬ 
erentially  dissolve,  we  plotted  the  final  masses  and  volumes  and  obtained 
densities  of  1.55  g  cm-3  for  IMX101, 1.67  g  cm-3  for  IMX104,  and  1.32  g  cnr 
3  for  PAX21. 
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Finally,  to  check  for  mass  closure,  we  compared  the  dissolved  masses 
measured  by  HPLC  with  the  difference  between  the  initial  and  ending 
masses  measured  by  analytical  balance  (Table  3).  For  IMX101  and  104,  we 
had  good  mass  balance  with  relative  errors  generally  less  than  5%.  The  er¬ 
rors  were  larger  for  PAX21  but  below  10%. 


Table  3.  Mass  balance  for  IMX101,  IMX104,  and  PAX21  particles  used  in  the  pCT  drip  dissolution  tests.  Also 

tabulated  are  the  calculated  volumes  for  each  particle.* 


Mass  (mg)  Scale 

Mass  (mg) 

Relative  Error 

Volume  (mm3) 

O 

II 

t=  End 

Difference 

HPLC 

(%) 

O 

II 

t=  End 

Difference 

IMX101 

#1 

18.97 

12.51 

6.46 

6.30 

2.48 

12.24 

8.23 

4.0 

#2 

18.65 

13.71 

4.94 

4.89 

1.01 

11.53 

8.58 

2.9 

#3 

24.03 

15.88 

8.15 

7.82 

4.05 

14.10 

10.09 

4.0 

#4 

20.37 

14.50 

5.87 

5.54 

5.62 

12.50 

9.61 

2.9 

A 

35.11 

19.91 

B 

16.58 

9.48 

C 

3.57 

1.97 

D 

36.81 

22.50 

p  (mg/mm3) 

1.67 

1.55 

IMX104 

#1 

18.18 

8.86 

9.32 

8.93 

4.17 

10.99 

5.47 

5.5 

#2 

22.37 

10.53 

11.84 

11.69 

1.30 

13.33 

6.03 

7.3 

#3 

52.57 

33.3 

19.27 

20.03 

3.95 

31.08 

19.52 

11.6 

#4 

37.32 

19.81 

17.51 

17.63 

0.67 

22.87 

11.94 

10.9 

p  (mg/mm3) 

1.67 

1.67 

PAX21 

#1 

33.02 

22.29 

10.73 

9.70 

9.60 

20.25 

17.64 

2.6 

#2 

16 

10.38 

5.62 

5.25 

6.58 

10.14 

7.7 4 

2.4 

#3 

15.77 

10.95 

4.82 

4.52 

6.22 

10 

7.99 

2.0 

#4 

30.75 

22.24 

8.51 

8.40 

1.29 

20.21 

16.24 

3.97 

B 

7.95 

4.98 

C 

82.08 

52.6 

p  (mg/mm3) 

1.51 

1.32 
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Figure  11.  Measurements  made  possible  by  pCT:  (a)  volume  of  particles  and  crystal 
constituents  in  IM  particles  and  (b)  Density  calculated  from  the  mass  and  pCT  volume 

measurements. 
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Volume  (mm3) 


2. 2. 2. 2  Dissolution  drip  tests 

We  used  well-controlled  laboratory  drip  experiments  to  measure  mass  dis¬ 
solved  as  a  function  of  water  volume  for  detonated  particles  of  IMXioi  and 
IMX104.  These  particles  were  dissolved  completely. 
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Figure  12  shows  results  for  the  drip  tests  on  four  IMX101  particles  (mass  of 
A  >  B  >  C  >  D).  To  show  the  compositional  variation  among  particles,  we 
plotted  the  ratio  of  each  compound  measured  to  its  initial  mass  (estimated 
based  on  the  mass  of  the  particle  and  the  putative  bulk  composition)  ver¬ 
sus  the  cumulative  water  volume.  The  smallest  particle  D,  for  example,  has 
more  NQ  and  less  NTO  and  DNAN  than  IMX101  bulk  composition,  where¬ 
as  particle  C  has  more  NTO  and  less  NQ  than  the  bulk.  The  sum  of  the  dis¬ 
solved  components  was,  nevertheless,  close  to  the  mass  of  the  original  par¬ 
ticle  (Table  4). 


Figure  12.  Percent  of  the  expected  mass  dissolved  for  compounds  in  IMX101  versus 

water  volume  (mL). 
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Table  4.  Mass  balance  for  IMX101  and  IMX104  drip  dissolution  tests.  The  expected  mass  was 
calculated  based  on  the  initial  mass  of  the  particle  and  the  putative  compositions  of  the 

formulations. 


IMX101 

IMX104 

NTO 

(mg) 

NQ 

(mg) 

DNAN 

(mg) 

Total 

(mg) 

NTO 

(mg) 

DNAN 

(mg) 

RDX 

(mg) 

Total 

(mg) 

A 

1 

expected 

13.58 

25.13 

29.21 

67.92 

5.16 

3.12 

1.46 

9.74 

measured 

12.26 

20.69 

26.8 

59.76 

5.55 

2.95 

1.67 

10.17 

%  error 

12 

4 

B 

2 

expected 

8.03 

14.86 

17.26 

40.15 

15.07 

9.10 

4.27 

28.44 

measured 

7.79 

13.94 

14.90 

36.63 

15.50 

8.32 

4.42 

28.24 

%  error 

9 

1 

c 

3 

expected 

5.74 

10.63 

12.35 

28.72 

13.84 

8.36 

3.92 

26.11 

measured 

6.56 

9.65 

12.39 

28.60 

14.18 

7.84 

4.07 

26.09 

%  error 

0 

0 

D 

4 

expected 

3.29 

6.08 

7.06 

16.43 

13.95 

8.42 

3.95 

26.32 

measured 

2.75 

7.20 

5.80 

15.75 

13.73 

7.93 

4.05 

27.18 

%  error 

4 

3 

These  plots  show  that  NTO  is  the  first  phase  dissolved,  followed  by  NQ, 
and  that  DNAN  is  last  component  dissolved,  consistent  with  the  results 
described  above  and  with  the  solubility  of  each  component  (Table  l).  The 
shape  of  the  mass  loss  curves  is  also  informative.  NTO  is  lost  quickly  as 
evident  by  a  steep  rise  in  mass  loss  followed  by  a  flat  plateau.  NQ  has  a  less 
steep  initial  rise  but  also  reaches  a  plateau,  although  for  the  largest  parti¬ 
cle,  A,  this  occurs  near  the  end  of  the  test.  The  DNAN  mass  loss  curves  are 
quasi  linear  for  A  and  D  and  are  described  by  two  lines  of  different  slope 
for  particles  B  and  C.  For  particles  B,  C,  and  D,  the  steepest  section  of  the 
DNAN  mass  loss  curve  occurs  before  NQ  has  completely  dissolved.  The 
dissolution  of  DNAN  tracks  that  of  NQ  in  this  region,  suggesting  that  NQ 
is  controlling  its  dissolution.  After  NQ  is  completely  dissolved,  DNAN  dis¬ 
solves  at  a  slower  rate.  We  think  that  NTO  dissolves  independently  and 
that  the  dissolution  of  NTO  and  NQ  leave  holes  in  the  particle  that  in¬ 
crease  DNAN’s  surface  area  and  its  dissolution.  As  soon  as  NQ  has  dis¬ 
solved,  the  DNAN  dissolves  independently. 
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The  mass  balances  for  these  tests  are  good  for  particles  C  and  D  and  are 
not  as  good  for  the  two  larger  particles,  A  and  B  (Table  4).  We  suspect  par¬ 
ticle  B  was  lost  close  to  the  end  of  the  test  because  the  funnel  holding  par¬ 
ticle  B  overflowed  after  1249  mL  of  effluent  had  been  collected  (day  115) 
after  which  the  concentration  for  DNAN  was  below  detection.  Although 
only  a  small  volume  of  water  overflowed,  we  suspect  the  particle  was  float¬ 
ing  on  the  surface  (these  IM  particles  can  trap  air  in  their  porous  struc¬ 
tures)  and  was  lost  along  with  the  water.  For  particle  A,  however,  the  fun¬ 
nel  did  not  overflow;  and  we  can  see  that  the  concentration  of  DNAN  in 
the  water  samples  decreased  gradually  until  it  was  below  the  detection 
limit.  In  this  case,  it  is  possible  that  a  small  piece  broke  from  the  particle 
between  weighing  the  sample  and  placing  it  in  the  frit  or  that  the  particle 
contained  small  pieces  of  metal  from  the  detonation  that  added  to  its  mass 
but  not  to  its  dissolved  mass. 

Figure  13  shows  results  for  the  drip  tests  on  four  IMX104  particles  (mass 
of  2  >  4  >  3  >  1).  Here  again,  NTO  dissolved  rapidly  and  independently  of 
the  DNAN-RDX  matrix.  The  percent  mass  loss  curves  for  the  DNAN- 
RDX  matrix  are  fairly  linear  with  a  larger  percentage  of  DNAN  dissolving 
early  in  the  test  and  a  larger  fraction  of  RDX  and  HMX  dissolving  later  in 
the  test.  Here  we  think  that  the  RDX  crystals  are  controlling  the  late  stage 
dissolution  of  the  particle  and  that  the  RDX  needs  to  either  dissolve  or  fall 
off  to  allow  the  underlying  DNAN  to  dissolve.  This  process  is  similar  to 
what  occurs  in  Comp  B  (Lever  et  al.  2005). 

Table  4  shows  the  mass  balances  for  IMX104  tests.  In  all  cases,  the  relative 
error  between  the  initial  mass  of  the  particle  and  the  mass  recovered  from 
the  water  samples  was  below  5%.  These  laboratory  drip  tests  indicate  that 
the  formulations  are  not  being  significantly  photo-  or  bio-transformed  be¬ 
cause  most  of  the  mass  (100  ±  5%)  is  recovered.  If  these  IM  formulations 
were  transforming  into  new  compounds  the  mass  balance  would  be  lower. 

Although  PAX21  is  no  longer  being  manufactured  and  used  in  munitions, 
for  completeness,  we  dripped  water  on  one  PAX21  particle.  Figure  14 
shows  the  mass  loss  as  a  function  of  water  volume.  The  pCT  data  do  not 
show  interconnections  between  the  AP  crystals  and  the  surface,  yet  the  da¬ 
ta  show  that  all  of  the  AP  dissolved  in  the  first  12  mL  of  water.  This  finding 
indicates  that  water  was  able  to  access  the  interior  of  this  51  mg  particle. 
The  data  clearly  show  that  the  AP  in  PAX21  dissolves  quickly  and  inde¬ 
pendently  of  the  DNAN-RDX  matrix.  The  PAX21  matrix  has  more  RDX, 
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with  a  2:1  RDX:DNAN  ratio,  than  does  the  IMX104  matrix,  with  a  1:1 
RDX:DNAN  ratio.  Consequently,  PAX21  particles  may  take  longer  to  dis¬ 
solve  than  comparable  IMX104  particles. 

Figure  13.  Percent  of  the  expected  mass  dissolved  for  compounds  in  IMX104  versus 
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Figure  14.  Drip  dissolution  test  of  a  51  mg  PAX21  particle. 


2.2.3  Outdoor  dissolution  tests 

We  set  up  the  outdoor  dissolution  tests  on  10  June  2013,  and  they  are  still 
ongoing.  We  have  collected  water  samples  on  eight  separate  occasions  and 
analyzed  six  sets  of  water  samples  (Figure  15).  A  sum  of  the  dissolved 
compound  masses  shows  that  the  IMX101  chunks  have  lost  between  21% 
to  65%  of  their  masses,  the  IMX104  chunks  between  33%  and  80%,  and 
the  PAX21  between  45%  and  71%  of  their  masses.  We  also  find  that  these 
formulations  are  much  more  friable  than  chunks  of  high  energetics.  In  less 
than  6  months,  8  of  the  12  IM  chunks  have  split;  and  all  of  the  chunks  have 
shed  millimeter-sized  particles  (Figure  16a,  b,  and  c).  The  average  volume 
of  rainfall  interacting  with  these  chunks  was  739  ±  82  mL.  For  compari¬ 
son,  of  the  12  chunks  of  TNT,  12  chunks  of  Comp  B,  5  chunks  of  tritonal, 
and  5  chunks  of  C4  (Composition  4,  a  common  variety  of  the  plastic  explo¬ 
sive  known  as  “C”)  placed  outdoors  to  weather  and  dissolve,  4  of  the 
chunks  split,  14  spalled  millimeter-sized  particles,  and  9  developed  cracks 
(Taylor  et  al.  2010)  in  three  years. 
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Figure  15.  Percent  mass  dissolved  of  the  mass  expected  for 
the  compounds  in  5  chunks  of  IMX101,  5  chunks  of 
IMX104,  and  2  chunks  of  PAX21  set  outside  to 
weather  and  dissolve. 
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Figure  16.  Sequential  images  of  IMX101  particle  #1  and  IMX104  particle  #3. 
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Our  tests  suggest  that  phototransformation  of  their  surfaces  is  occurring 
because  (1)  the  IM  pieces  changed  color  from  cream  or  white  to  orange  or 
brick  red,  (2)  an  unknown  peak  develops  in  the  HPLC  where  polar  sub¬ 
stances  elute,  (3)  this  peak  is  largest  in  samples  collected  during  the  sum¬ 
mer  months  when  sun  irradiance  is  highest,  and  (4)  indoor  (no  sunlight) 
tests  on  IM  explosives  showed  good  mass  balance  (100  ±  5%)  and  no  un¬ 
accounted  for  HPLC  peaks.  Although  we  have  not  yet  calculated  a  mass 
balance  for  these  samples  (this  will  be  done  at  the  end  of  the  test)  we  ex¬ 
pect  that  it  will  be  low  because  of  the  transformation  of  the  IM  compounds 
into  unknown  products  that  we  are  not  analyzing  for. 

2.2.4  pH  of  the  IM  solutions 

NTO  is  known  to  be  acidic  in  solution  (Chipen  et  al.  1966;  Le  Campion  et 
al.  1997),  and  we  measured  the  pH  of  the  both  the  drip  and  the  outdoor 
water  samples  to  quantify  how  this  changed  during  dissolution  of  the  par¬ 
ticle.  For  the  IMX101  drip  samples,  we  measured  pH  values  initially  in  the 
low  3  range  that  increased  to  neutral  as  the  NTO  concentration  decreased 
to  10  mg  L1  (Figure  17).  The  pH  of  the  IMX104  samples  showed  a  similar 
trend  with  the  difference  being  that,  because  IMX104  contains  propor¬ 
tionally  more  NTO  than  does  IMX101,  the  concentrations  were  higher  and 
the  pH  levels  lower  for  the  first  few  IMX104  samples.  Figure  13  shows  that 
there  is  a  good  correlation  between  NTO  concentration  and  pH  values  be¬ 
tween  2  and  5  for  all  samples.  At  NTO  concentrations  below  20  mg  L_1,  the 
pH  can  vary  between  4  and  6;  and  at  NTO  concentrations  below  10  mg  L1, 
the  pH  of  the  solution  is  close  to  neutral.  The  color  of  the  solution  was  also 
found  to  correlate  with  the  NTO  concentration— the  most  yellow  colored 
solutions  had  the  highest  NTO  concentrations  and  the  lowest  pH  values. 
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Figure  17.  Plot  showing  the  relationship  between  NTO  concentration 
and  pH  for  indoor  drip  samples  and  outdoor  samples.  These  data 
look  similar  to  a  titration  curve. 
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3  Phototransformation  studies 


3.1  Methods 


3.1.1  Phototransformation  of  NTO  and  DNAN  in  solution 

To  evaluate  direct  phototransformation  rates  in  aqueous  solutions  of  NTO 
and  DNAN,  we  followed  ASTM  Method  E  896-92  ( Standard  Test  Method 
for  Conducting  Aqueous  Direct  Photolysis  Tests )  (ASTM  International 
2005).  We  exposed  solutions  of  NTO  and  DNAN  to  the  simulated  sunlight 
and  quantified  the  concentration  changes  using  an  HPLC  with  UV/Vis  de¬ 
tector.  Specifically,  we  dissolved  a  sample  in  deionized  water,  adjusted  to 
pH  7  using  a  phosphate  buffer,  to  achieve  a  1  mg  L_1  concentration  and 
placed  it  in  quartz  NMR  5  mm  tubes.  We  conducted  phototransformation 
measurements  by  using  a  solar  simulator,  SUNTEST  CPS+  Xenon  Test  In¬ 
strument  by  Atlas  Material  Testing  Technology  (Chicago,  IL)  (Figure  18). 
The  Atlas  Suntest  CPS+  uses  an  air-cooled  xenon  lamp  to  simulate  UV  and 
visible  solar  radiation  in  the  300-800  nm  range.  Irradiance  (the  radiative 
flux,  or  the  amount  of  light  incident  on  a  unit  area  of  a  surface)  of  765 
W/m2  was  used  for  a  period  of  200-500  min  with  samples  collected  and 
analyzed  over  time. 


Figure  18.  SUNTEST  CPS+  Xenon  Test  Instrument  used  to 
evaluate  phototransformation  of  NTO  and  DNAN  in  solution. 


We  measured  concentrations  of  DNAN  and  NTO  in  solution  before  and 
after  irradiation  by  using  an  Agilent  1200  series  HPLC  with  C18  column  at 
30°C,  a  10  pL  manual  injection  loop,  and  acetonitrile/ water  (60/40)  mo¬ 
bile  phase.  A  UV/Vis  detector  measured  DNAN  at  300  nm  (using  360  nm 
as  a  reference)  and  NTO  at  375  nm  (using  390  nm  as  a  reference).  No 
analysis  of  phototransformation  products  was  done  during 
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phototransformation  rate  studies.  Dark  samples  of  DNAN  and  NTO  were 
monitored  as  the  control.  We  found  both  compounds  to  be  stable  in  the 
dark  over  the  time  of  the  experiment. 

According  to  ASTM  896  (ASTM  International  2005),  phototransformation 
of  ionizable  materials  can  exhibit  marked  pH  effects  that  are  attributable 
to  changes  in  speciation.  Therefore,  if  different  species  are  present  in  the 
pH  5-9  range,  testing  should  be  conducted  at  2-3  pH  values  throughout 
the  pH  5-9  range.  NTO  is  an  ionizable  compound  (Figure  19)  with  a  pKa  of 
3.67  (Kofman  et  al.  1980).  A  pH  of  1  ppm  NTO  solution  in  deionized  water 
equals  to  a  pH  of  4.98.  Therefore,  in  this  study,  direct 
phototransformation  of  NTO  was  measured  at  pH  1.5,  4.2,  7.0,  and  8.8. 

We  achieved  these  pH  values  by  using  phosphate,  acetate,  phosphate,  and 
borate  buffers,  respectively. 

Figure  19.  Ionized  and  neutral  forms  of  NTO  (Smith  and  Cliff  1999). 


1  N 

In  addition  to  direct  phototransformation,  we  measured  the  indirect 
phototransformation  of  neat  DNAN  and  NTO  in  solution.  To  determine 
the  contribution  of  dissolved  natural  organic  matter  present  in  solution  to 
the  transformation  of  IM  compounds  under  light,  we  used  U.S.  EPA  meth¬ 
od  OPPTS  835.5270,  “Indirect  Photolysis  Screening  Test”  (U.S.  EPA  1998). 
Humic  acid  (from  Sigma  Aldrich,  St.  Louis,  MO)  in  60  ppm  dissolved  OC 
concentration  served  as  natural  organic  matter. 

We  took  initial  phototransformation  rate  measurements  (including  those 
at  varying  pH  for  NTO)  at  35°C  because  it  is  the  lowest  temperature  that 
the  solar  simulator  can  maintain.  This  temperature  is  within  the  recom¬ 
mended  standard  range  of  test  temperatures  for  phototransformation  ac¬ 
cording  to  ASTM  E896-92  (ASTM  International  2005)  but  is  higher  than 
can  be  expected  in  some  environments.  Consequently,  to  describe  temper¬ 
ature  dependence  of  phototransformation  rates,  we  performed  direct 
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3.2 


phototransformation  measurements  at  pH  7  for  four  different  tempera¬ 
tures  between  70°C  and  35°C.  Control  solutions  set  in  the  dark  indicated 
that  both  compounds  were  stable  across  the  temperature  range  used.  Once 
we  obtained  rates  at  different  temperatures,  we  used  the  following  equa¬ 
tion  to  relate  measured  reaction  rates  to  the  activation  energy  of  the  reac¬ 
tion: 


_Ea 

kphoto  ^6  RT  (1) 

where  reaction  rate,  kphoto,  is  a  function  of  temperature,  T,  and  activation 
energy  of  the  reaction,  Ea. 

Calculated  activation  energy  of  the  phototransformation  of  each  com¬ 
pound  allows  us  to  determine  the  rate  reaction  at  any  temperature.  The 
HYDRUS  model  is  able  to  incorporate  information  about  the  energy  of  ac¬ 
tivation  (Ea)  for  the  phototransformation  reactions  for  DNAN  and  NTO 
into  predictions  of  IM  fate  in  the  environment.  We  did  all  measurements 
in  triplicate  and  observed  them  to  be  highly  reproducible  with  error  gener¬ 
ally  smaller  than  the  size  of  the  symbols  on  the  plot. 

Results 

3.2.1  Direct  and  indirect  phototransformation  of  NTO 

Figure  20  presents  NTO  phototransformation  as  a  function  of  time.  Meas¬ 
ured  NTO  concentrations  were  highly  reproducible  with  standard  devia¬ 
tion  of  the  mean  less  than  1%.  All  phototransformation  reactions  exhibited 
strong  first-order  kinetics  (Figure  20).  The  transformation  rate  constants 
were  estimated  by  fitting  regression  lines  to  plots  of  the  logarithm  of  NTO 
concentration  as  a  function  of  time.  The  regression  lines  have  high  correla¬ 
tion  coefficients  (R2= 0.99),  indicating  a  good  fit  and  high  confidence  in  the 
estimated  rate  constants. 

The  pH  of  the  solution  strongly  affected  the  phototransformation  of  NTO, 
with  rates  being  highest  at  extreme  pH  values  and  lowest  at  neutral  pH 
(Figure  20a;  Table  5).  The  rate  constant  at  pH  7  was  estimated  to  be  equal 
to  0.0014  min-1  or  0.084  hr-1,  resulting  in  a  half-life  of  495  min  or  8.3  hr. 
Measured  NTO  phototransformation  rate  constants  were  higher  than  con¬ 
stants  measured  previously  for  RDX  (Bordeleau  et  al.  2013),  which  NTO  is 
replacing  in  explosive  formulations. 
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Figure  20.  Phototransformation  of  NTO  as  affected  by  (a)  solution  pH  and 
(b)  the  presence  of  natural  organic  matter. 
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NTO  phototransformation  was  strongly  affected  by  the  presence  of  natural 
organic  matter  (humic  acids)  with  the  phototransformation  rate  being  al¬ 
most  three  times  higher  in  the  presence  of  organic  matter  compared  to 
rates  measured  when  no  organic  matter  was  present  (Figure  20b,  Table  5). 

3.2.2  Direct  and  indirect  phototransformation  of  DNAN 

DNAN  also  experienced  first-order  kinetics  of  phototransformation 
(Figure  21).  At  pH  7,  the  phototransformation  rate  constant  was  0.0021 
min-1,  resulting  in  a  half-life  of  330  min.  For  comparison,  the  half-life  of 
TNT  in  pure  water  is  700  min  and  200  min  in  fresh  water.  The  shorter 
half-life  of  fresh  water  may  be  due  to  the  presence  of  natural  organic  mat¬ 
ter  (O’Sullivan  et  al.  2011;  Prak  and  O’Sullivan  2012).  Unlike  TNT  and 
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NTO,  which  have  been  shown  to  be  affected  by  organic  matter  in  solution, 
DNAN’s  phototransformation  was  not  influenced  by  the  presence  of  or¬ 
ganic  matter,  with  a  phototransformation  rate  of  0.0020  min-1  and  a  half 
life  of  346  min,  values  that  were  not  significantly  different  from  ones 
measured  for  DNAN  in  water  without  organics  present. 

Figure  21.  Phototransformation  rate  of  DNAN  in  water  and  in  water  with 
humic  acid.  Error  bars  are  smaller  than  the  symbols. 
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Measured  DNAN  phototransformation  rate  constants  were  similar  to  ones 
determined  in  Rao  et  al.  (2013),  which  reported  an  average  half-life  for 
DNAN  of  0.70  days  with  values  ranging  between  0.11  and  1.51  days.  The 
range  in  values  depended  on  the  wavelength  and  the  intensity  of  the  light 
source,  the  initial  DNAN  concentration,  and  the  temperature.  For  condi¬ 
tions  similar  to  ones  used  in  our  study,  Rao  et  al.  (2013)  estimated  a  0.25 
day  half  life. 


Table  5.  Determined  phototransformation  rate  constants,  kphoto,  and  half-lives,  fe,  for  NTO 
and  DNAN  measured  at  35°C,  765  Watts/m2. 


NTO 

DNAN 

pH 

kphoto 

(min-1) 

kphoto 

(h-1) 

& 

(min) 

& 

(h) 

tu 

(day) 

kphoto 

(min-1) 

kphoto 

(h-1) 

tu 

(min) 

tu 

(h) 

tu 

(day) 

Direct 

1.5 

0.0049 

0.294 

139 

2.3 

0.10 

4.2 

0.0026 

0.156 

267 

4.5 

0.19 

7.0 

0.0014 

0.084 

495 

8.3 

0.34 

0.0021 

0.126 

330 

5.5 

0.23 

8.8 

0.0051 

0.306 

138 

2.3 

0.10 

Indirect 

7.0 

0.0041 

0.246 

168 

2.8 

0.12 

0.0020 

0.120 

346 

5.8 

0.24 
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3.2.3  Phototransformation  of  NTO  and  DNAN  as  a  function  of  temperature 

In  addition,  we  performed  a  set  of  experiments  to  determine  the  effect  of 
temperature  on  phototransformation  rate  constants  of  NTO  and  DNAN 
and  to  calculate  the  Ea  of  phototransformation  reactions  for  these  com¬ 
pounds.  We  found  that  phototransformation  of  NTO  was  relatively  insen¬ 
sitive  to  temperature  (Figures  22,  23,  and  24;  Table  6),  with  tx/2  equal  to 
0.30  days  at  35°C  and  0.27  days  at  70°C,  resulting  in  an  Ea  of  2.65  kJ 
mol-1  and  a  predicted  half-life  of  0.31  days  at  27°C.  Low  Ea  for  NTO  sug¬ 
gests  catalysis  by  decomposition  products. 


Figure  22.  Effect  of  temperature  on  DNAN  and  NTO  phototransformation  rates. 
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Figure  23.  Natural  logarithm  of  the  phototransformation  rate  constant  plotted  as 
a  function  of  temperature  for  DNAN  and  NTO. 
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Phototransformation  of  DNAN  was  temperature  dependent  (Figure  22,  23, 
24;  Table  6)  with  Ea  =  27.8  kJ  mol-1;  DNAN  had  a  measured  ti/2  of  0.13 
days  at  70°C,  0.23  days  at  35°C,  and  a  predicted  half-life  of  0.50  days  at 
27°C.  Determined  activation  energies  for  DNAN  were  lower  than  deter¬ 
mined  by  Rao  et  al.  (2013),  which  reported  an  Ea  of  53  kJ  mol-1.  Knowing 
the  activation  energies  of  the  phototransformation  reactions  allows  us  to 
calculate  the  transformation  rates  for  different  environments. 
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Figure  24.  Measured  and  predicted  values  of  NTO  and  DNAN  half-lives  as  a 
function  of  temperature. 
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Table  6.  Phototransformation  rate  constant,  kphoto,  and  half-life,  fe,  of  NTO  and  DNAN 
tabulated  as  a  function  of  temperature.  Measurements  were  made  at  765  W/m2  irradiance. 


NTO 

DNAN 

kphoto 

*y2 

kphoto 

*y2 

*y2 

Temperature 

(min-1) 

(min) 

(h) 

(day) 

(min-1) 

(min) 

(h) 

(day) 

-X 

o 

o 

CN 

0.00155 

446 

7.4 

0.31 

0.00091 

763 

12.7 

0.53 

00 

CJl 

o 

o 

0.00159 

437 

7.3 

0.30 

0.00120 

576 

9.6 

0.40 

00 

o 

o 

0.00164 

424 

7.1 

0.29 

CJl 

o 

o 

o 

0.00170 

407 

6.8 

0.28 

0.00208 

333 

5.6 

0.23 

O) 

o 
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0.00174 

398 

6.6 

0.28 

0.00271 

256 

4.3 

0.18 

->J 

o 

o 
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0.00176 

395 

6.6 

0.27 

0.00361 

192 

3.2 

0.13 

^Extrapolated. 
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4  Soil  interactions 

4.1  Materials  and  methods 

4.1.1  Soil  sample  collection  and  characterization 

To  quantify  IM-soil  interactions,  we  used  soils  having  a  wide  range  of 
properties  and  that  were  collected  on  military  installations  located  across 
the  United  States.  The  U.S.  Army  National  Guard  Environmental  Direc¬ 
torate  provided  uncontaminated  soil  samples  from  a  number  of  different 
locations.  In  addition,  we  used  soils  that  were  previously  collected  studied 
for  their  interaction  with  explosive  and  propellant  constituents  (Dontsova 
et  al.  2006,  2007,  2009a,  2009b;  Taylor  et  al.  2012).  Using  the  same  soils 
allows  us  to  compare  our  results  for  IMs  to  those  obtained  for  other  ener¬ 
getics.  We  used  only  surface  soils  (top  30  cm). 

We  air  dried,  sieved  (<  2  mm),  and  characterized  the  collected  soils  prior 
to  the  experiments.  We  analyzed  OC  by  dry  combustion  according  to  the 
U.S.  EPA’s  standard  operating  procedure  for  sediments  (U.S.  EPA  2010). 
To  ensure  a  complete  removal  of  inorganic  carbon,  we  treated  a  known 
mass  of  soil  with  5%  phosphoric  acid  (H3P04).  Following  acid  pretreat¬ 
ment,  we  used  a  Shimadzu  SSM-5000A  Solid  Sample  Combustion  Unit  to 
analyze  OC.  A  Beckman  Coulter  LS  13  320  Laser  Diffraction  Particle  Size 
Analyzer  (by  Beckman  Coulter,  Inc.,  Fullerton,  CA)  analyzed  particle  sizes. 
A  Beckman  Coulter  SA  3100  measured  specific  surface  area  (SSA)  by  N2 
adsorption,  applying  Brunauer-Emmett-Teller  theory  (Brunauer  et  al. 
1938).  We  used  a  VWR  sympHony  SB70P  Benchtop  Digital  pH  and  EC 
(electrical  conductivity)  Meter  to  test  a  1:1  soiksolution  mixture. 

4.1.2  IM  compounds  and  their  transformation  products 

For  the  experiments,  we  obtained  the  technical  grade  NTO  and  DNAN 
from  the  U.S.  Army  Armament  Research,  Development  and  Engineering 
Center  (ARDEC),  Picatinny  Arsenal.  Calibration  standards  for  NTO  were 
prepared  using  the  same  NTO,  while  DNAN  standards  were  prepared  us¬ 
ing  DNAN  (CAS:  119-27-7,  98%  purity)  purchased  from  Alfa-Aesar  (Ward 
Hill,  MA).  Known  daughter  compounds  of  DNAN  (Perreault  et  al.  2011; 
Olivares  et  al.  2013),  2-methoxy-5-nitroaniline  (MENA)  (CAS:  99-59-2, 
98%  purity)  and  2,4-Diaminoanisole  (CAS:  615-05-4  analytical  standard), 
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that  we  used  to  calibrate  HPLC  were  obtained  from  Sigma- Aldrich  (St. 
Louis,  MO). 

4.1.3  Analytical  methods 

We  quantified  NTO  concentration  by  using  a  Dionex  Ultimate  3000  HPLC 
equipped  with  a  diode  array  detector  (from  ThermoFisher,  MA).  We 
adapted  the  running  method  from  Le  Campion  et  al.  (1999a).  The  mobile 
phase,  made  from  AcN  with  0.1%  TFA  and  deionized  water  (75:25),  was 
run  isocractically  at  1  mL  min-1.  Oven  temperature  was  set  at  32°C.  NTO 
and  its  transformation  products  were  separated  using  a  Thermo  Scientific 
Hypercarb  Column.  We  detected  NTO  at  approximately  2.5  min  by  using  a 
315  nm  wavelength.  The  UV  detector  was  set  at  220  nm  to  monitor  for 
presence  of  ATO  (Le  Campion  et  al.  1998),  but  we  did  not  detect  ATO  in 
any  analyzed  samples. 

The  parameters  used  for  analyzing  DNAN  samples  were  the  following:  a 
mobile  phase  ratio  of  43:57  methanol  and  water  and  flow  rate  of  1  mL 
min-1.  The  column  used  for  analysis  was  the  Thermo  Scientific  Acclaim 
reversed  phase  C-18  column.  The  wavelengths  used  were  300  nm  for 
DNAN,  254  nm  for  MENA,  and  210  nm  for  2,4-diaminoanisole  (Olivares 
et  al.  2013). 

4.1.4  Batch  soil  adsorption  studies 

To  determine  the  fate  of  IM  compounds  in  soils  and  the  effect  of  soil  het¬ 
erogeneity  on  soil  interactions  with  IM  compounds,  we  conducted  a  series 
of  kinetic  and  equilibrium  batch  soil  sorption  experiments  (Roy  et  al. 
1992).  In  these  experiments,  we  mixed  soils  with  solutions  of  NTO  or 
DNAN  and  allowed  them  to  equilibrate  for  a  predetermined  amount  of 
time.  We  then  separated  the  soil  from  the  solution  and  analyzed  the  su¬ 
pernatant  to  determine  the  amount  of  NTO  or  DNAN  remaining  in  solu¬ 
tion.  We  conducted  preliminary  experiments  to  determine  the  optimal 
soiksolution  ratios  of  NTO  and  DNAN  (Table  7)  and  the  effect  of  back¬ 
ground  salt  solution  on  adsorption  of  these  compounds.  For  follow-on 
studies,  we  selected  tests  that  resulted  in  7096-90%  of  the  original  NTO  or 
DNAN  amount  remaining  in  the  solution.  The  experiments  used  deionized 
water  or  0.005  mol  L-1  CaCl2  solution  to  establish  the  effect  of  background 
salt  solution  on  adsorption. 
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To  create  a  l  mg  L1  input  solution,  we  added  10  mL  of  a  too  mg  L1  stock 
solution  of  NTO  to  990  mL  of  deionized  water  or  of  0.005  mol  L_1  CaCl2 
solution.  We  diluted  DNAN  further  to  0.5  mg  L-1.  Soils  were  weighed  into 
clean,  labeled  30  mL  borosilicate  glass  centrifuge  tubes  according  to  Table 
7.  We  then  pipetted  20  mL  of  1  mg  L_1  NTO  (or  0.5  mg  L_1  DNAN)  in  water 
or  0.005M  CaCl2  solution  into  each  centrifuge  tube.  An  aliquot  of  the  input 
solution  was  collected  and  stored  to  determine  the  initial  concentration  of 
the  solute  later.  We  conducted  duplicate  or  triplicate  tests,  depending  on 
observed  variability,  and  prepared  blanks  by  putting  the  input  solution 
(prepared  either  in  water  or  0.005M  CaCl2)  into  a  centrifuge  tube  holding 
no  adsorbent.  A  reciprocating  shaker  agitated  blanks  and  samples  contin¬ 
uously  at  110  rpm  for  24  ±  0.5  hr.  Experiments  were  conducted  under  con¬ 
stant  conditions,  away  from  light  and  at  room  temperature  (22  ±  3°C). 

After  24  hr  of  agitation,  we  recorded  any  changes  observed  in  the  adsor¬ 
bent  or  solutions  and  centrifuged  the  tubes  for  20  min  at  4700  rpm  (4816 
relative  centrifugal  force).  A  minimum  of  3  mL  of  supernatant  was  filtered 
through  0.45  pm  Millex-HV  PVDF  filter  (from  EMD  Millipore,  Darmstadt, 
Germany)  and  placed  in  labeled  4  mL  amber  vials.  For  experiments  with 
NTO,  we  pipetted  0.1  mL  of  this  sample  into  labeled  autosampler  vials  for 
analysis  and  diluted  them  to  1  mL  with  75:25  AcN:  deionized  water  mix¬ 
ture.  We  analyzed  the  supernatant  solution  for  target  compounds,  NTO 
and  NTO  transformation  product,  ATO,  by  HPLC,  and  analyzed  solutions 
containing  DNAN  by  using  HPLC  directly  after  filtering  without  additional 
dilution.  Section  4.1.3  of  this  report  describes  HPLC  methods  for  NTO, 
DNAN,  and  their  transformation  products.  The  pH  and  EC  of  each  sample 
was  determined  on  an  aliquot  using  a  VWR  SympHony  SB70P  Benchtop 
Digital  pH  and  EC  Meter  (VWR,  Radnor,  PA).  We  wrapped  all  prepared 
solutions  in  aluminum  foil  (to  prevent  phototransformation)  and  sealed 
them  with  parafilm,  storing  samples  and  stock  solutions  in  a  refrigerator 
at  approximately  4°C. 

We  observed  no  significant  difference  in  NTO  or  DNAN  soil  interaction 
from  samples  prepared  with  or  without  a  background  solution  of  0.005 
mol  L_1  CaCl2.  All  remaining  tests  used  a  0.005  mol  L_1  CaCl2  background 
solution.  We  selected  the  soiksolution  ratio  indicating  between  10%  and 
30%  adsorption  and  used  it  in  kinetic  adsorption  and  transformation  stud- 
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Table  7.  Soihsolution  ratios  tested  in  NTO  and  DNAN  soil  adsorption  studies. 


Soil:Solution  Ratio 
(g  ml.-1) 

Soil  Weight 
(g) 

Volume  of  Solution  Containing  Solute 
(mL) 

1:4 

5 

20 

1:10 

2 

20 

1:20 

1 

20 

1:40 

0.5 

20 

1:60 

0.33 

20 

The  optimal  equilibrium  soihsolution  ratio  for  NTO  in  all  soils  was  1:4 
while  for  DNAN  1:60  ratio  worked  best  for  majority  of  the  soils.  To  observe 
changes  in  NTO  and  DNAN  soil  adsorption  and  transformation  with  time, 
we  conducted  a  kinetic  adsorption  study.  Here  we  prepared  in  triplicate 
treatment  samples  and  blanks  without  soil  and  equilibrated  them  for  1,  4, 
8,  24,  48,  72,  96,  and  120  hr;  we  then  centrifuged,  filtered,  and  stored 
them  as  previously  described.  We  tested  aliquots  for  pH  and  EC  and  de¬ 
termined  the  time  interval  for  attaining  chemical  equilibrium  for  each  soil 
type.  Equilibrium  time  was  defined  as  the  minimum  time  required  to 
reach  a  rate  of  change  in  solute  concentration  of  less  than  5%  in  24  hr  (Roy 
et  al.  1992).  NTO  and  DNAN  transformation  rate  constants,  k,  were  de¬ 
termined  from  the  change  in  solute  concentration  over  time,  dC/dt.  For 
first-order  rate,  the  equation  takes  the  following  form 

-dC/dt=kC.  (2) 

The  integrated  form  of  the  first-order  rate  law  is 

lnC=-kt+lnCo  (3) 

where  Co  is  initial  concentration  of  target  analyte  in  solution.  The  slope  of 
this  plot  was  used  to  determine  k. 

After  determining  the  equilibrium  time,  we  constructed  soil  adsorption 
isotherms.  Isotherms  comprised  eight  input  concentrations:  0.78, 1.56, 
3.13,  6.25, 12.5,  25,  50,  and  100  mg  Lr1  NTO  or  2.03,  4.06,  8.13, 16.25, 

32.5,  65,  and  130  mg  Lr1  DNAN  in  0.005  mol  Lr1  CaCl2. 

We  determined  adsorption  coefficients  from  measured  solution  concentra¬ 
tions  and  calculated  sediment  concentrations  using  linear  and  Freundlich 
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adsorption  isotherms.  The  following  equation  describes  the  linear  iso¬ 
therm: 


5=  KdC 


(4) 


where  S  is  the  contaminant  concentration  in  soil  (mg  kg-1),  Kd  is  the  ad¬ 
sorption  coefficient  (mL  g_1),  and  C  is  the  contaminant  concentration  in 
solution  (mg  Lr1).  The  Freundlich  isotherm  is  described  by 


5=  Kf&  (5) 

where  Kf  is  the  Freundlich  adsorption  coefficient  and  n  is  an  empirical  pa¬ 
rameter  that  indicates  the  affinity  of  the  compound  for  the  adsorbent.  We 
used  the  regression  analysis  function  in  Windows  Excel  2010  to  determine 
linear  distribution  coefficients  (Kd),  confidence  intervals  for  the  estimate, 
R2,  and  probabilities  of  parameters  being  significantly  different  from  zero. 
The  Freundlich  equation  was  used  in  linearized  form  to  facilitate  applica¬ 
tion  of  linear  regression  tools  to  determine  adsorption  parameters.  Among 
the  soils,  we  compared  parameters  from  the  equation  (linear  of 
Freundlich)  that  best  described  the  data  (largest  R2)  for  the  majority  of  the 
soils  to  elucidate  the  effect  of  soils  on  adsorption.  We  related  determined 
parameters  to  measured  soil  properties,  such  as  their  particle  size  distribu¬ 
tion,  organic  carbon  content,  and  pH,  to  determine  the  effect  of  geochemi¬ 
cal  heterogeneity  on  adsorption  of  IM  compounds  and  to  help  predict  the 
environmental  fate  of  these  compounds  in  soils. 

4.1.5  NTO  transport  studies 

To  measure  the  transport  of  NTO  in  soils,  we  conducted  saturated  column 
transport  studies  under  steady-state  and  transient  conditions  (Dontsova  et 
al.  2006,  2009b).  Solutions  of  NTO  were  used  as  the  source.  We  selected 
some  of  the  soils  used  in  the  batch  studies  to  use  in  column  transport  ex¬ 
periments.  In  the  column  effluent,  we  determined  concentrations  of  NTO 
(presence  of  ATO  was  monitored  but  not  observed)  and  a  conservative 
tracer  and  used  these  to  construct  breakthrough  curves. 

We  conducted  saturated  flow  column  experiments  by  using  Supelco  (Belle- 
fonte,  PA)  glass  tubes  (7  cm  in  length  with  a  1.18  cm  internal  diameter) 
with  PTFE  (polytetrafluoroethylene)  caps.  Small  columns  minimize  the 
amount  of  source  material  needed  and  also  reduce  hazardous  waste  pro¬ 
duced  in  the  study.  We  packed  11-12  g  of  soil  homogenously  over  a  bottom 
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layer  of  silanized  glass  wool  that  prevented  migration  of  particles.  A  layer 
of  glass  wool  overlay  the  packed  soil  and  acted  to  evenly  distribute  incom¬ 
ing  flow.  To  avoid  air  entrapment,  we  saturated  columns  from  the  bottom 
up  with  a  0.005  mol  L1  CaCl2  solution.  After  saturation,  we  determined 
pore  volume  as  the  volume  of  solution  required  to  saturate  the  packed  col¬ 
umn.  We  attached  Tygon  microbore  tubing,  running  from  a  Cole-Parmer 
(Vernon  Hills,  IL)  Master  flex  peristaltic  pump,  to  the  top  of  each  column 
to  supply  NTO  or  DNAN  solution  in  CaBr2  background  at  a  given  flux.  We 
used  a  0.02  mL  min-1  flow  rate  or  1.1  cm  hr-1  flux.  We  continuously  col¬ 
lected  outflow  (effluent)  into  4  mL  amber  vials  by  using  a  Teledyne  ISCO 
(Lincoln,  NE)  Foxy  200  Fraction  collector  with  a  200  vial  capacity.  A  max¬ 
imum  of  ten  columns  could  run  simultaneously.  After  about  6  to  16  pore 
volumes,  depending  on  the  soil,  inflow  was  switched  to  0.005  mol  L1 
CaCl2  solution;  and  the  flow  continued  for  another  6-16  pore  volumes  to 
evaluate  NTO  desorption  from  the  soil.  We  analyzed  a  conservative  tracer 
(bromide,  Br)  by  using  Ion  Chromatography  (Dionex  ICS  5000  with  diode 
array)  to  distinguish  between  non-equilibrium  attributed  to  physical  ver¬ 
sus  chemical  processes  and  to  determine  dispersion  and  diffusion  in  the 
columns. 

As  the  same  soils  were  used  in  batch  adsorption  and  transformation  stud¬ 
ies,  we  were  able  to  compare  determined  parameters,  allowing  us  to  verify 
the  processes  responsible  for  NTO  attenuation.  In  addition,  as  soils  from 
some  of  these  sites  have  been  used  in  previous  explosive  and  propellants 
studies  (Dontsova  et  al.  2006,  2007 , 2008,  2009a,  2009b;  Taylor  et  al. 
2012),  we  were  able  to  compare  the  transport  of  IM  compounds  to  those  of 
other  energetics. 

We  performed  flow-interruption  tests,  wherein  flow  was  stopped  for  24  hr, 
for  some  experiments  to  quantify  mass-transfer  processes,  such  as  diffu¬ 
sion  and  sorption  (Brusseau  et  al.  1989, 1997).  Stopping  the  flow  allows 
the  system  to  equilibrate.  This  is  reflected  in  decreases  in  concentration  if 
degradation  or  kinetic  sorption  is  taking  place. 

4.1.6  IM  formulation  dissolution  and  transport  studies 

To  study  dissolution  and  transport  of  IMs,  we  used  undetonated  particles 
of  IMX101,  IM  formulation  containing  DNAN  and  NTO,  and  detonation 
residues  as  the  source  (Dontsova  et  al.,  2006,  2009b).  Camp  Guernsey  and 
Camp  Swift  soils  were  used  in  column  transport  and  IM  solid  dissolution 
experiments.  For  the  experiments  employing  the  formulations,  we  placed 
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explosives  particles  on  the  surface  of  soil  saturated  with  0.005  mol  L1 
CaCl2  and  started  the  flow  of  the  0.005  mol  Lr1  CaBr2  solution  (with  Br 
serving  as  the  conservative  tracer).  Once  steady-state  outflow  was  reached, 
we  measured  the  desorption  of  explosives  from  the  soil  by  removing  the 
explosives  particles  and  continuing  the  flow  of  the  0.005  M  CaCl2  solution 
until  we  no  longer  detected  the  target  analytes  in  the  outflow. 

We  collected  and  analyzed  outflow  from  the  columns,  the  water  that  inter¬ 
acted  with  the  explosives  and  then  with  soil,  by  using  a  Dionex  UHPLC  (ul¬ 
tra  high  performance  liquid  chromatography)  with  Diode  Array  Detector 
to  evaluate  dissolution,  degradation,  and  adsorption  of  IM  compounds  and 
their  transformation  products.  The  Br  tracer  was  used  to  distinguish  be¬ 
tween  physical  and  chemical  non-equilibrium  processes  and  to  determine 
dispersion  and  diffusion  in  the  columns.  We  measured  concentration  of 
the  Br-  tracer  by  using  ion  chromatography  (Dionex  ICS  5000  ion  chro¬ 
matograph  with  conductivity  detector,  ThermoFisher,  CA). 

Collecting  effluent  directly  into  the  vials  (Figure  25)  prevented  contact  be¬ 
tween  the  target  solute  and  the  plastic  tubing,  eliminating  possible  adsorp¬ 
tion  pathways  that  can  have  significant  impact  on  measured  solution  con¬ 
centrations. 

We  determined  concentrations  of  NTO,  DNAN,  their  transformation 
products,  and  NQ,  which  is  also  present  in  IMX101,  in  the  column  effluent 
and  used  these  to  construct  breakthrough  curves. 

We  used  a  0.01  mL  min-1  flow  rate  (0.55  cm  hr-1  flux)  to  match  drip  stud¬ 
ies.  No  flow  interruption  was  done  during  the  IMX101  dissolution  experi¬ 
ments. 
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Figure  25.  Experimental  setup  for  the  IM  column  transport  studies. 


4.2  Results 

4.2.1  Soil  sample  collection  and  characterization 

We  used  eleven  different  soils  for  this  study.  Ten  of  these  soils  were  col¬ 
lected  at  military  training  ranges,  and  Catlin  silt  loam  was  collected  on  the 
University  of  Illinois  at  Urbana-Champaign  university  farm. 

The  Catlin  and  two  additional  soils,  Plymouth  sandy  loam  and  Sassafras 
loam,  were  previously  used  to  study  soil  interactions  and  transport  of  pro¬ 
pellant  constituents;  studies  of  high  explosives  had  also  used  Catlin  and 
Plymouth  soils.  New  soils  were  sent  to  us  by  the  Army  National  Guard- 
Environmental  Directorate  from  training  sites  affiliated  with  the  U.S.  Ar¬ 
my  National  Guard  (Table  8).  Facilities  sampled  included  AZ  National 
Guard  Florence  Military  Reservation;  Camp  Butner,  Army  National  Guard 
Training  Center,  NC;  Camp  Swift  Army  National  Guard  Training  Center, 
TX;  Camp  Guernsey,  a  Premier  Joint  Training  Center  in  Guernsey,  WY, 
serving  both  a  U.S.  Air  Force  and  National  Guard.  Additionally,  samples 
came  from  two  locations  in  Montana:  Fort  William  Henry  Harrison,  the 
Montana  National  Guard  training  facility,  where  small  arms,  heavy  ma¬ 
chine  gun,  and  pyrotechnics  are  fired,  and  Limestone  Hills  Training  Area, 
where  artillery  and  grenade  training  and  open  burn/open  detonation  is 
done.  These  soils  belong  to  the  same  series  but  contain  different  amount  of 
organic  carbon,  twice  as  much  for  Fort  Harrison  compared  to  Limestone 
Hills.  Samples  were  also  collected  on  the  small  arms  range  at  Arnold  Air 
Force  Base,  Tullahoma,  TN. 
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Soils  that  were  collected  earlier  are  identified  hereafter  using  their  series 
name  while  new  soils  are  identified  using  the  location  where  they  were  col¬ 
lected.  Table  8  lists  the  soil  series  name  and  classification  for  all  soils. 


Table  8.  Soils  used  to  study  interaction  of  IMs  with  the  soils.  Series  names  and  classifications 
are  according  to  U.S.  Taxonomy  and  location  where  collected. 


Soil  name 

Soil  classification 

Location 

Catlin 

Catlin  silt  loam,  mixed,  mesic,  superactive 

Oxyaquic  Argiudoll 

Urbana,  IL 

Fort  Harrison 

Musselshell  sandy  loam,  carbonatic,  frigid  Aridic 
Calciustepts 

Fort  William  Henry 
Harrison,  MT 

Arnold  AFB 

Captina  silt  loam,  siliceous,  mesic,  typic 

Fragiudults 

Arnold  Air  Force  Base 
(AFB),  TN 

Plymouth 

Plymouth  loamy  sand,  mesic,  coated  Typic 
Quartzipsamment 

Massachusetts  Military 
Reservation,  MA 

Camp  Butner 

Helena  sandy  loam,  mixed,  semiactive,  thermic 
Aquic  Hapludults 

Camp  Butner,  NC 

Limestone  Hills 

Musselshell  sandy  loam,  carbonatic,  frigid  Aridic 
Calciustepts 

Limestone  Hills,  MT 

Sassafras 

Sassafras  loam,  siliceous,  mesic  Typic  Hapudult 

Aberdeen  Proving 

Ground,  MD 

Camp  Gruber 

Verdigris  silt  loam,  mixed,  superactive,  thermic 
Cumulic  Hapludolls 

Camp  Gruber,  OK 

Camp  Guernsey 

Keeline-Turnercrest  loam,  mixed,  superactive, 
calcareous,  mesic  Ustic  Torriorthents 

Camp  Guernsey,  WY 

Florence  MR 

Cherioni  loam,  mixed,  superactive,  hyperthermic, 
shallow  Typic  Haplodurids 

Florence  Military 
Reservation  (MR),  AZ 

Camp  Swift 

Bergstrom  sandy  clay  loam,  mixed,  superactive, 
thermic  Cumulic  Haplustolls 

Camp  Swift,  TX 

Selected  soils  represent  a  range  of  climate  and  conditions  on  military 
ranges  where  live-fire  training  is  occurring  (Figure  26).  They  encompass 
multiple  soil  orders,  including  Entisols,  Inceptisols,  Mollisols,  Aridisols, 
and  Ultisols,  with  different  moisture  regimes,  from  very  dry  to  wet.  Table  9 
summarizes  the  physical  and  chemical  properties  for  the  11  soils  selected 
for  this  study.  All  studied  soils  had  mixed  particle  size,  ranging  from  fine 
(clay  loam)  to  coarse  (loamy  sand). 
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Figure  26.  Map  of  the  continental  Unites  States  indicating  location  of  soil  collection. 
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Table  9.  Measured  physical  and  chemical  properties  of  soils  used  in  adsorption  and  transport 
studies  with  NTO,  DNAN,  and  IM  formulations. 


Soil 

Texture 

Clay 

(%) 

Silt 

(%) 

Sand 

(%) 

pH* 

ECt 

(pS  cm-1) 

SSA* 

(msg-1) 

oc§ 

(%) 

Catlin 

silt  loam 

25.6 

65.5 

8.9 

7.31 

492 

6.4 

5.28 

Fort  Harrison 

sandy  loam 

8.7 

36.5 

54.9 

6.67 

449 

7.4 

3.88 

Arnold  AFB 

silt  loam 

11.4 

65.5 

23.1 

6.66 

131 

7.8 

2.68 

Plymouth 

loamy  sand 

4.4 

20.4 

75.2 

4.23 

206 

1.7 

2.45 

Camp  Butner 

sandy  loam 

7.7 

25.9 

66.4 

6.69 

219 

4.8 

2.42 

Limestone  Hills 

sandy  loam 

11.2 

35.7 

53.1 

7.54 

539 

10.5 

1.99 

Sassafras 

loam 

16.4 

42.3 

41.4 

4.40 

212 

7.17 

1.30 

Camp  Gruber 

clay  loam 

32.3 

44.9 

22.8 

5.39 

74 

38.3 

0.83 

Camp  Guernsey 

loam 

4.1 

12.5 

83.4 

8.21 

477 

3.9 

0.77 

Florence  MR 

loam 

26.8 

33.5 

39.7 

8.00 

417 

33.0 

0.45 

Camp  Swift 

sandy  clay 
loam 

23.7 

20.8 

55.6 

7.83 

203 

15.1 

0.34 

*ln  l:lsoil:water 
t  EC  =  Electrical  conductivity 
*  SSA  =  specific  surface  area 
§  OC  =  organic  carbon 
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4.2.2  NTO  batch  adsorption  studies 

We  designed  our  preliminary  experiments  to  establish  the  effect  of  a  back¬ 
ground  electrolyte,  0.005  mol  L1  CaCl2  solution,  on  sorption  of  NTO  in 
soils  and  to  determine  the  soiksolution  ratio  that  would  result  in  10%- 
30%  retention  of  NTO  on  the  soil.  An  adsorption  of  io%-30%  results  in 
the  smallest  relative  error  in  adsorption  experiments  (Roy  et  al.  1992; 
McDonald  and  Evangelou  1997),  and  a  1:4  soiksolution  ratio  came  closest 
to  the  target  solute  removal  for  NTO.  NTO  is  negatively  charged  under  the 
environmental  pH  range;  and  therefore,  electrolytes  can  influence  NTO 
sorption.  However,  preliminary  experiments  indicated  that  there  was  no 
effect  of  background  electrolytes  on  sorption  at  the  concentration  used. 

Kinetic  experiments  indicated  that  the  majority  of  soils  reached  equilibri¬ 
um  with  NTO  in  solution  at  24  hr  with  the  exception  of  Catlin  and  Plym¬ 
outh  soils,  which  reached  equilibrium  at  8  hr,  and  Limestone  Hills  soil, 
which  required  48  hr.  Kinetic  experiments  also  allowed  us  to  determine 
transformation  rate  constants  for  NTO  in  soils  (Figure  27;  Table  10).  Con¬ 
trol  NTO  solution  (blank)  did  not  indicate  a  decrease  in  concentration  over 
120  hr  kinetic  sorption  experiments.  When  soils  were  present,  the  majority 
of  soils  had  a  faster  decrease  in  solution  concentrations  over  the  first  24  hr 
followed  by  a  slower  decrease.  We  detected  no  organic  products  of  NTO 
transformation.  Figure  27  presents  results  for  the  kinetic  experiment  using 
Plymouth  soil,  including  a  plot  of  concentrations  normalized  to  input  over 
time  and  the  natural  logarithm  of  this  value  plotted  as  a  function  of  time. 
The  slope  of  the  linear  regression  between  time  and  natural  logarithm  of 
the  concentration  is  used  to  determine  the  rate  of  the  first-order  transfor¬ 
mation  reaction.  Estimates  ranged  between  0.0004  hr-1  in  Camp  Guern¬ 
sey  and  0.0221  hr-1  in  Catlin.  Rate  estimates  were  highly  significant,  99%, 
and  different  from  zero  for  all  soils,  except  for  Limestone  Hills  where  it 
was  only  significant  at  the  95%  level. 

Measured  transformation  rates  were  lower  for  soils  with  less  OC,  such  as 
Camp  Guernsey,  Florence  MR,  and  Camp  Swift.  There  was  a  positive  rela¬ 
tionship  between  the  percent  OC  present  in  the  soil  and  the  measured 
transformation  rate  constant  k  (probability,  P  =  0.02),  however,  R2  was 
low  at  0.46  (Figure  28).  Other  soil  parameters,  such  as  clay  content,  pH, 
and  SSA,  poorly  correlated  with  k.  The  measured  NTO  transformation  rate 
in  Plymouth  soil,  0.0043  hr-1,  was  about  three  times  slower  than  one  ob¬ 
served  previously  for  RDX  (0.013  hr-1)  (Dontsova  et  al.  2006). 
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Figure  27.  Plot  of  (a)  NTO  solution  concentrations  normalized  by  input  over 
time  in  the  presence  of  Plymouth  soil  and  (b)  the  same  concentrations  plotted 
as  the  natural  logarithm  of  concentrations  in  order  to  determine  the 
transformation  rate  constant  of  NTO.  Error  bars  equal  one  standard  error  of 
the  mean  (n  =  9  with  3  independent  replicate  test  tubes  and  3  repeat  HPLC 
measurements  for  each  sample). 
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NTO  adsorbed  very  weakly  to  the  studied  soils  as  can  be  expected  for  a 
negatively  charged  compound  in  a  matrix  that  also  possesses  a  net  nega¬ 
tive  charge.  All  determined  adsorption  coefficient  {Kci)  values  were  less 
than  l  cm3  g-i  with  the  largest  observed  Kd  equal  to  0.51  cm3  g-i  in  Camp 
Gruber  soil.  For  most  soils,  both  Freundlich  and  linear  adsorption  iso¬ 
therm  described  the  observed  adsorption  well;  but  Freundlich  usually  re¬ 
sulted  in  a  slightly  better  fit  (Figures  29-32;  Table  10).  The  Freundlich  pa¬ 
rameter  n  was  not  similar  among  the  soils,  with  the  majority  smaller  than 
one,  but  was  above  one  for  the  two  soils  with  the  highest  OC  content,  Cat- 
lin  and  Fort  Harrison.  Fits  to  both  linear  and  Freundlich  isotherms  were 
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worse  for  soils  that  had  very  little  adsorption  (Camp  Guernsey,  Florence 
MR,  and  Camp  Swift  had  Kd  less  than  o.i  cm3  g-i  and  Camp  Butner  0.12 
cm3  g-i)  due  to  a  larger  percentage  of  error  for  these  measurements. 

Table  10.  Fate  and  transport  parameters  for  NTO  in  studied  soils:  Freundlich  adsorption 
parameters,  /frand  rr,  soil  adsorption  coefficient,  Kd,  adsorption  coefficient  normalized  to 
fraction  of  organic  carbon  in  soils,  Kod,  transformation  rate  constant,  k,  and  Z?2  values  for  the 

fits. 


Soil 

Kf 

n 

R2 

Kd 

(cm3 

g-1) 

R2 

Koo 

(cm3  g-1) 

k 

(h-i) 

R2 

Catlin 

0.21 

1.03 

0.94 

0.21 

0.92 

3.98 

0.0221 

0.84 

Fort  Harrison 

0.27 

1.07 

0.98 

0.35 

0.95 

9.02 

0.0021 

0.98 

Arnold  AFB 

0.58 

0.86 

0.98 

0.34 

0.94 

12.69 

0.0044 

0.94 

Plymouth 

0.82 

0.89 

0.99 

0.50 

0.96 

20.41 

0.0043 

0.97 

Camp  Butner 

0.74 

0.54 

0.77 

0.12 

0.72 

4.96 

0.0021 

0.98 

Limestone  Hills 

0.33 

0.88 

0.92 

0.21 

0.92 

10.55 

0.0123 

0.60 

Sassafras 

0.90 

0.86 

0.99 

0.48 

0.96 

36.92 

0.0080 

0.99 

Camp  Gruber 

0.54 

0.99 

1.00 

0.51 

0.99 

61.45 

0.0025 

0.98 

Camp  Guernsey 

0.06 

0.48 

0.18 

0.02 

0.21 

2.60 

0.0004 

0.75 

Florence  MR 

0.09 

0.77 

0.55 

0.06 

0.57 

13.33 

0.0005 

0.91 

Camp  Swift 

0.10 

0.84 

0.92 

0.04 

0.59 

11.76 

0.0009 

0.93 

Figure  28.  Correlation  between  measured  NTO  transformation  rate  constants 
(A)  and  soil  0C.  P  =  0.02250. 
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Figure  29.  Adsorption  isotherms  calculated  for  NTO  in  Catlin,  Fort 
Harrison,  and  Limestone  Hills  soils.  The  red  dashed  line  indicates  the 
linear  adsorption  isotherm  fit  to  the  measured  adsorption  data 
(equation  is  in  red)  while  the  solid  black  line  is  a  fit  of  the  Freundlich 
isotherm  with  the  equation  presented  in  black. 
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Figure  30.  Adsorption  isotherms  calculated  for  NTO  in  Sassafras, 
Plymouth,  and  Camp  Gruber  soils.  The  red  dashed  line  indicates  the 
linear  adsorption  isotherm  fit  to  the  measured  adsorption  data 
(equation  is  in  red)  while  the  solid  black  line  is  a  fit  of  the  Freundlich 
isotherm  with  the  equation  presented  in  black. 
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Figure  31.  Adsorption  isotherms  calculated  for  NTO  in  Camp  Butner, 
Florence  MR,  and  Camp  Guernsey  soils.  The  red  dashed  line  indicates 
the  linear  adsorption  isotherm  fit  to  the  measured  adsorption  data 
(equation  is  in  red)  while  the  solid  black  line  is  a  fit  of  the  Freundlich 
isotherm  with  the  equation  presented  in  black. 
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Figure  32.  Adsorption  isotherms  calculated  for  NTO  in  Arnold  AFB  and 
Camp  Swift  soils.  The  red  dashed  line  indicates  the  linear  adsorption 
isotherm  fit  to  the  measured  adsorption  data  (equation  is  in  red)  while 
the  solid  black  line  is  a  fit  of  the  Freundlich  isotherm  with  the  equation 

presented  in  black. 


0  20  40  00  80  100  120 


0  20  40  60  80  100  120 

C  (mg  L-1} 


Soil  pH  was  the  strongest  indicator  of  NTO  affinity  for  soil  surfaces  (Figure 
33).  There  was  a  highly  significant  negative  relationship  between  linear 
adsorption  coefficients  and  soil  pH.  We  observed  no  relationship  between 
Kd  values  and  OM,  clay,  or  SSA.  The  K/- values  had  similar  if  slightly  weaker 
( R 2  =  0.7818,  P  =  0.00030)  negative  linear  relationships  with  pH. 
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Figure  33.  Measured  NTO  adsorption  coefficients  ( Kd )  plotted  against 
soil  pH.  P  =  0.00011. 


For  Catlin  soil  that  has  a  high  pH  and  a  high  OC  content,  NTO  was  ab¬ 
sorbed  about  10  times  less  than  RDX,  0.21  versus  2.03  cm3  g-i  (Dontsova 
et  al.  2009a),  while  for  Plymouth  soil  (low  pH,  less  OC)  NTO  was  more 
strongly  adsorbed  than  RDX,  0.82  versus  0.65  cm3  g-i  (Dontsova  et  al. 
2006). 

4.2.3  DNAN  batch  adsorption  studies 

We  found  that  the  best  soilrsolution  ratio  for  the  majority  of  DNAN  ad¬ 
sorption  experiments  was  1:60.  We  observed  no  effect  of  the  background 
electrolyte  solution  (0.005  mol  L1  CaCl2)  on  adsorption.  Kinetic 
adsorption  studies  conducted  to  determine  the  best  equilibration  time  for 
obtaining  isotherms  and  to  determine  the  transformation  rate  of  DNAN  in 
soils  indicated  that  the  concentration  of  DNAN  in  the  control  (blank— no 
soil)  samples  did  not  change  over  120  hr  (Figure  34).  For  samples  with 
soil,  however,  there  was  a  sharp  drop  in  DNAN  concentrations  in  the  first 
hours  followed  by  a  steady  decrease  in  concentration.  This  change 
followed  first-order  kinetics. 

Transformation  rate  constants,  k,  ranged  between  0.0006  and  0.0073  hr-1 
for  Camp  Swift  and  Camp  Gruber  soils,  resectively.  Although  k  appears  to 
increase  with  OC,  the  trend  is  not  statistically  significant  (P  =  0.78101). 
Our  study  found  that  other  soil  parameters,  such  as  pH,  clay  content,  and 
SSA,  are  not  correlated  with  k.  No  known  products  of  DNAN 
tranformation  were  measured  in  the  solutions.  All  soils  reached  equilibri¬ 
um  at  24  hr,  except  Sassafras  loam,  which  reached  equilibrium  with 
DNAN  at  48  hr. 
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Figure  34.  Plot  of  (a)  normalized  DNAN  solution  concentrations  in 
contact  with  Camp  Guernsey  soil  and  (b)  the  same  concentrations 
plotted  as  the  natural  logarithm  whose  slope  gives  the  transformation 

rate,  /f,  of  NTO. 
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DNAN  transformation  was  significantly  slower  than  observed  for  TNT.  For 
Plymouth  soil,  the  TNT  transformation  rate  constant  was  0.21  hr-1 
(Dontsova  et  al.  2006)  while  the  DNAN  constant  was  0.0070  hr-1. 

For  all  studied  soils,  the  Freundlich  isotherm  resulted  in  a  highly  signifi¬ 
cant  fit  to  the  data.  Isotherms  were  similar  in  shape  (Figures  35-38).  Pa¬ 
rameter  n,  an  exponent  in  the  Freundlich  equation,  was  close  for  different 
soils,  estimates  averaging  0.76  with  a  range  between  0.56  and  0.97.  How¬ 
ever,  95%  confidence  intervals  of  n  estimates  overlapped  for  the  majority 
of  the  soils,  indicating  that  they  were  not  significantly  different  from  each 
other  at  95%.  For  three  soils,  Sassafras,  Camp  Guernsey  and  Camp  Swift,  n 
was  not  significantly  different  from  1  (resulting  in  a  linear  isotherm). 
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Figure  35.  Adsorption  isotherms  calculated  for  DNAN  in  Arnold  AFB,  Camp 
Swift,  and  Sassafras  soils.  The  red  dashed  line  indicates  the  linear 
adsorption  isotherm  fit  to  the  measured  adsorption  data  (equation  is  in  red) 
while  the  solid  black  line  is  a  fit  of  the  Freundlich  isotherm  with  the 
equation  presented  in  black. 
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Figure  36.  Adsorption  isotherms  calculated  for  DNAN  in  Catlin,  Florence 
MR,  and  Camp  Gruber  soils.  The  red  dashed  line  indicates  the  linear 
adsorption  isotherm  fit  to  the  measured  adsorption  data  (equation  is  in 
red)  while  the  solid  black  line  is  a  fit  of  the  Freundlich  isotherm  with  the 
equation  presented  in  black. 
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Figure  37.  Adsorption  isotherms  calculated  for  DNAN  in  Camp  Guernsey, 
Camp  Butner,  and  Plymouth  soils.  The  red  dashed  line  indicates  the 
linear  adsorption  isotherm  fit  to  the  measured  adsorption  data  (equation 
is  in  red)  while  the  solid  black  line  is  a  fit  of  the  Freundlich  isotherm  with 
the  equation  presented  in  black. 
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Figure  38.  Adsorption  isotherms  calculated  for  DNAN  in  Limestone  Hills 
and  Fort  Harrison  soils.  The  red  dashed  line  indicates  the  linear 
adsorption  isotherm  fit  to  the  measured  adsorption  data  (equation  is  in 
red)  while  the  solid  black  line  is  a  fit  of  the  Freundlich  isotherm  with  the 
equation  presented  in  black. 
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Linear  regressions  were  also  highly  significant  with  slopes  {Kci  values)  sig¬ 
nificantly  different  from  zero.  While,  for  the  majority  of  soils,  a  Freundlich 
isotherm  described  DNAN  adsorption  the  best,  for  four  soils,  Fort  Harri¬ 
son,  Sassafras,  Camp  Guernsey,  and  Camp  Swift,  a  linear  isotherm  result¬ 
ed  in  a  better  fit  to  the  data  (Table  li).  Normalizing  estimated  Kd  values  to 
OC  content  in  the  soil  ( Koc )  resulted  in  a  decrease  in  the  percent  standard 
deviation  of  the  estimates  (from  6i%  to  49%). 
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Table  11.  Fate  and  transport  parameters  for  DNAN  in  studies  soils:  Freundlich  adsorption 
parameters,  Kf  and  nr,  linear  adsorption  coefficient,  Kd,  adsorption  coefficient  normalized  to 
fraction  of  organic  carbon  in  soils,  Kod,  transformation  rate  constant,  k,  and  /^values  for  the 

fits. 


Soils 

Kr 

n 

R2 

Kd 

(cm3  g-1) 

R2 

Koo 

(cm3  g-1) 

k 

(h-i) 

R2 

Catlin 

34.00 

0.62 

0.98 

5.95 

0.92 

112.7 

0.0047 

0.88 

Fort  Harrison 

13.86 

0.77 

0.9 

6.32 

0.93 

162.9 

0.0010 

0.45 

Arnold  AFB 

14.4 

0.68 

0.93 

3.39 

0.78 

126.5 

0.0022 

0.92 

Plymouth 

10.05 

0.83 

0.98 

4.38 

0.94 

178.6 

0.0070 

0.51 

Camp  Butner 

15.35 

0.56 

0.91 

2.05 

0.92 

84.8 

0.0018 

0.67 

Limestone  Hills 

10.75 

0.77 

0.96 

4.96 

0.92 

250 

0.0019 

0.78 

Sassafras 

2.26 

0.97 

0.85 

1.89 

0.72 

228.3 

0.0013 

0.54 

Camp  Gruber 

7.62 

0.72 

0.84 

1.99 

0.91 

240.6 

0.0073 

0.29 

Camp  Guernsey 

1.72 

0.85 

0.93 

0.93 

0.97 

121.4 

0.0041 

0.85 

Florence  MR 

6.59 

0.74 

0.98 

1.91 

0.94 

424.4 

0.0007 

0.35 

Camp  Swift 

1.27 

0.84 

0.7 

0.6 

0.84 

175.5 

0.0006 

0.96 

Determined  DNAN  adsorption  coefficients,  both  Kd  and  Kf,  resulted  in 
highly  significant  positive  relationships  with  percent  OC  in  the  soil  (Figure 
39).  No  other  measured  soil  properties,  such  as  clay  content,  pH,  or  SSA, 
correlated  with  Kt  or  Kf  values. 

Comparing  DNAN  with  TNT  adsorption  coefficients  for  Catlin  and  Plym¬ 
outh  soils  showed  higher  adsorption  of  DNAN  in  Plymouth  soil  (Ki  of  4.38 
cm3  g-i  for  DNAN  and  0.63  to  1.6  cm3  g-i  for  TNT,  depending  on  the 
method  used)  (Dontsova  et  al.  2006)  but  less  adsorption  of  DNAN  in  Cat¬ 
lin  soil.  Although  the  Freundlich  parameters  were  similar  for  two  com¬ 
pounds  in  the  Catlin  soil  ( K/=  34.00,  n  =  0.62  for  DNAN  and  K/=  32.67,  n 
=  0.60  for  TNT),  the  Kd  values  were  higher  for  TNT  (17.9  vs.  5.95) 
(Dontsova  et  al.  2009a). 
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Figure  39.  Correlation  between  (a)  measured  DNAN  linear  (/O)  and  (b) 
Freundlich  (/G)  adsorption  coefficients  and  soil  percent  OC.  P  =  0.00061 
for  /fry  and  0.00011  for  Kf,  both  highly  significant  relationships. 
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4.2.4  NTO  transport  studies 

For  a  subset  of  the  eight  soils  used  in  the  batch  experiments,  we  conducted 
column  experiments  using  NTO.  All  experiments  were  conducted  in  dupli¬ 
cate  and  replicated  well.  Figures  40  through  43  show  breakthrough  curves 
for  Camp  Swift,  Camp  Guernsey,  Sassafras,  and  Arnold  AFB  soils.  The  ma¬ 
jority  of  the  soils  behaved  as  expected  based  on  the  batch  studies.  For  ex¬ 
ample,  the  early  breakthrough  curves  and  steady-state  concentrations  seen 
in  Camp  Swift  soil  (Figure  40)  were  consistent  with  low  Kct  and  k  estimated 
in  batch  studies  (Table  10),  while  Sassafras  breakthrough  curves  (Figure 
42)  were  consistent  with  higher  Kd  and  k  measured  in  batch  studies. 
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Figure  40.  Breakthrough  curves  for  NTO  in  Camp  Swift  soil.  The  dashed 
grey  vertical  line  indicates  timing  of  24  hr  flow  interruption.  The  tan 
vertical  line  indicates  the  time  when  the  solution  was  changed  back  to  a 
saturating  solution  to  observe  the  desorption  phase  of  the  isotherm. 
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Figure  41.  Breakthrough  curves  for  NTO  in  Camp  Guernsey  soil.  The 
dashed  grey  vertical  line  indicates  timing  of  24  hr  flow  interruption.  The  tan 
vertical  line  indicates  the  time  when  the  solution  was  changed  back  to  a 
saturating  solution  to  observe  the  desorption  phase  of  the  isotherm. 
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Figure  42.  Breakthrough  curves  for  NTO  in  Sassafras  soil.  The  dashed  grey 
vertical  line  indicates  tinning  of  24  hr  flow  interruption.  The  tan  vertical  line 
indicates  the  time  when  the  solution  was  changed  back  to  a  saturating 
solution  to  observe  the  desorption  phase  of  the  isotherm. 
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However,  several  soils,  Arnold  AFB  (Figure  43),  Catlin,  and  Limestone 
Hills  had  breakthrough  curves  that  showed  increases  in  transformation 
rate  with  time  and  no  recovery  after  the  interrupted  flow.  Both  trends  in¬ 
dicate  that  the  mechanism  of  transformation  may  be  changing  as  the  col¬ 
umn  becomes  increasingly  anaerobic  with  time. 


Figure  43.  Breakthrough  curves  for  NTO  in  Arnold  AFB  soil.  The  dashed 
grey  vertical  line  indicates  timing  of  24  hr  flow  interruption.  The  tan 
vertical  line  indicates  the  time  when  the  solution  was  changed  back  to  a 
saturating  solution  to  observe  the  desorption  phase  of  the  isotherm. 
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Table  12  presents  transport  parameters  determined  from  batch  and  col¬ 
umn  experiments.  As  Table  12  shows,  a  first-order  transformation  rate 
constant,  k,  and  values  estimated  using  batch  and  column  experiments 
were  significantly  different  from  each  other  for  all  soils  that  were  used  in 
both  types  of  experiments.  Their  95%  confidence  intervals  did  not  overlap. 
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Table  12.  Estimates  of  fate  and  transport  parameters;  linear  adsorption  coefficients,  /0,and  transformation 
rate  constants,  k,  for  NTO  determined  from  batch  and  column  experiments. 

CE  =  95%  confidence  interval. 


Soil 

Batch  estimates 

Column  estimates 

Kd  (cm3  g"1) 

k(h-i) 

Kd  (cm3  g-1) 

k(h-i) 

Estimate 

Lower 

CE 

Upper 

CE 

Estimate 

Lower 

CE 

Upper 

CE 

Estimate 

Lower 

CE 

Upper 

CE 

Estimate 

Lower 

CE 

Upper 

CE 

Catlin 

0.21 

0.18 

0.24 

0.0221 

0.0113 

0.0330 

0.46 

0.37 

0.56 

NA 

NA 

NA 

Arnold  AFB 

0.34 

0.29 

0.38 

0.0044 

0.0032 

0.0055 

0.83 

0.36 

1.29 

NA 

NA 

NA 

Camp  Butner 

0.10 

0.07 

0.14 

0.0021 

0.0018 

0.0024 

0.42 

0.20 

0.64 

0.185 

0.114 

0.255 

Limestone 

Hills 

0.21 

0.17 

0.26 

0.0123 

0.0023 

0.0223 

0.21 

0.13 

0.30 

NA 

NA 

NA 

Sassafras 

0.46 

0.43 

0.50 

0.0080 

0.0076 

0.0085 

0.86 

0.47 

1.24 

0.107 

0.091 

0.123 

Camp 

Guernsey 

0.02 

0.01 

0.04 

0.0004 

0.0002 

0.0006 

0.17 

0.13 

0.20 

0.066 

0.060 

0.072 

Florence  MR 

0.05 

0.03 

0.08 

0.0005 

0.0004 

0.0007 

0.21 

0.13 

0.28 

0.069 

0.067 

0.072 

Camp  Swift 

0.04 

0.01 

0.06 

0.0009 

0.0007 

0.0012 

0.05 

-0.01 

0.12 

0.003 

-0.001 

0.007 

Correlation  between  batch  and  column  NTO  linear  adsorption  coefficients 
(Kd)  was  better  (Figure  44)  (P  =  0.00185,  or  highly  significant).  However, 
the  slope  of  the  regression  was  higher  than  1,  indicating  that  there  was  a 
trend  towards  higher  Kd  estimates  in  the  column  experiments.  Confidence 
interval  overlapped  for  only  two  soils,  Limestone  Hills  and  Camp  Swift, 
indicating  that  estimates  were  different  at  95%. 


Figure  44.  Correlation  between  batch  and  column  NTO  linear  adsorption  coefficients  (Kd). 
P  =  0.00185  for  Kd,  which  is  highly  significant. 
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4.2.5  IMX101  dissolution  and  transport 

We  performed  IMX101  dissolution  and  transport  studies  using  Camp  Swift 
and  Camp  Guernsey  soils  that  were  used  earlier  in  the  batch  experiments 
with  NTO  and  DNAN.  Both  soils  were  also  part  of  the  NTO  column  stud¬ 
ies,  and  we  selected  them  for  the  dissolution  and  transport  studies  because 
of  their  low  OM  content  as  DNAN  adsorption  was  expected  to  be  the  time- 
limiting  step  in  the  experiment.  Soils  most  significantly  differed  in  their 
particle  size.  Camp  Guernsey  is  a  sandy  soil  while  Camp  Swift  has  high 
clay  content  (Table  9). 

IMX101  contains  43.5%  DNAN,  36.8%  NQ,  and  19.7%  NTO.  Breakthrough 
curves  were  observed  for  all  IMX  components.  In  addition,  we  detected 
MENA,  a  DNAN  transformation  product,  in  the  column  outflow.  MENA 
concentrations  in  the  effluent  were  more  than  100  times  smaller  than 
measured  DNAN  concentrations.  We  did  not  observe  MENA  in  any  of  the 
previous  experiments  with  DNAN  conducted  in  this  study.  However,  the 
presence  of  MENA  is  consistent  with  the  much  higher  concentrations  of 
parent  compound,  DNAN,  observed  in  IMX101  column  effluents. 

Behavior  of  IM  constituents  NTO,  DNAN,  and  NQ  was  consistent  with 
previously  known  (for  NQ)  and  determined  (for  NTO  and  DNAN)  fate  and 
transport  parameters  and  agreed  with  patterns  of  IMX  dissolution  ob¬ 
served  in  drip  studies  (Figures  45  and  46).  NTO  was  the  first  compound  to 
break  through.  NTO  also  had  the  highest  observed  outflow  concentrations 
for  all  experiments  despite  representing  only  about  20%  of  the  total 
IMX101  mass.  High  observed  concentrations  agree  with  the  high  solubility 
of  NTO  and  with  observations  in  drip  studies.  Early  NTO  breakthrough  is 
consistent  with  low  measured  adsorption  coefficients  for  both  studied  soils 
(0.02  and  0.04  cm3  g-i  for  Camp  Guernsey  and  Camp  Swift,  respectively). 
Nitroguanidine,  similarly  to  NTO  is  a  relatively  inert  compound  in  soil  en¬ 
vironments.  Nitroguanidine  Kd  values  for  Sassafras,  Plymouth,  and  Catlin 
soils  were  0.60,  0.44,  and  0.24  cm3  g-i  (Taylor  et  al.,  2012)  compared  to 
0.48,  0.50,  and  0.21  cm3  g-i  measured  for  NTO  in  the  same  soils.  Break¬ 
through  curves,  however,  indicated  a  slightly  later  arrival  of  NQ  compared 
to  NTO  and  a  delay  in  the  desorption  phase,  indicating  that  NQ  has  a 
higher  affinity  for  studied  soils  than  does  NTO.  Peak  concentrations  of  NQ 
were  lower  than  those  of  NTO  despite  being  present  at  higher  concentra¬ 
tions  in  solid  IMX101.  The  solubility  of  NTO  (16,642  mg  L_1)  is  about  four 
times  higher  than  that  of  NQ  (3800  mg  L1)- 
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Figure  45.  Breakthrough  curves  for  NTO;  NQ;  DNAN;  and  DNAN 
transformation  product,  MENA,  from  IMX101  particle  dissolution  in 
Camp  Swift  soil:  (a)  undetonated  particles  and  (b)  residues  of  low  order 
detonation.  The  flow  rate  is  0.01  mL  min-1.  The  black  vertical  line  marks 
the  time  when  the  IMX101  particles  were  removed  from  the  soil  surface 
to  observe  desorption  phase  of  the  breakthrough. 
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Figure  46.  Breakthrough  curves  for  NTO;  NQ;  DNAN;  and  DNAN 
transformation  product,  MENA,  from  undetonated  IMX101  particle 
dissolution  in  Camp  Guernsey  soil.  Flow  rate  is  0.01  ml_  min-1.  The 
black  vertical  line  marks  the  time  when  the  IMX101  particle  was 
removed  from  the  soil  surface  to  observe  the  desorption  phase  of  the 

breakthrough. 

Camp  Guernsey,  undetonated  IMX  101,  column  7 
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DNAN  was  the  last  IM  constituent  to  breakthrough.  It  also  had  corre¬ 
sponding  delay  in  the  desorption  phase  of  the  experiment.  In  addition, 
DNAN  exhibited  the  lowest  solution  concentrations  among  IM  constitu¬ 
ents  present  in  IMX101.  MENA  appeared  in  the  outflow  at  about  the  same 
time  as  DNAN  for  Camp  Swift  soil.  However,  in  Camp  Guernsey  soil  it 
appeared  earlier.  If  we  assume  that  transformation  of  DNAN  to  MENA  is 
occuring  in  solution,  we  can  conclude  that  MENA  is  adsorbing  to  soils  less 
than  DNAN.  There  was  also  a  difference  between  adsorption  of  MENA  in 
Camp  Swift  and  Camp  Guernsey  soils.  MENA  concentrations  remained 
high  at  the  end  of  the  experiment  in  Camp  Swift  but  started  to  descrease  in 
Camp  Guernsey  soil.  MENA’s  lower  affinity  to  Camp  Guernsey  soil  as 
compared  to  DNAN  exhibited  itself  in  a  higher  sensitivity  of  MENA  to 
changes  in  IMX101  dissolution,  such  as  greater  roughness  of  the 
breakthrough  curve  and  sharper  peaks. 


We  added  approximately  the  same  amount  of  IMX101  (0.05  g)  to  each 
column.  IMX101  in  columns  3,  4,  7,  and  8  for  the  set  of  experiments  per¬ 
formed  at  0.01  mL  min-1  was  undetonated  while  particles  used  in  columns 
1,  2, 5,  and  6  came  from  low-order  detonations.  Differences  in  column  out¬ 
flow  between  undetonated  IMX  particles  and  residues  of  low-order 
IMX101  detonations  (Figure  45)  were  within  particle-to-particle  variability 
for  undetonated  IMX101. 
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5  Summary 

The  solubility  of  the  individual  constituents  of  IM  formulations  differ  by 
orders  of  magnitude,  producing  porous  particles  as  they  dissolve,  which 
increases  both  the  surface  areas  and  the  dissolution  rates  of  the  particles. 
The  drip  dissolution  tests  and  the  pCT  data  show  that  the  components  in 
these  explosive  formulations  dissolve  sequentially  and  in  the  order  pre¬ 
dicted  by  their  solubility.  The  good  mass  balances  measured  for  the  drip 
tests  indicate  that  the  dissolved  formulations  do  not  significantly  photo-  or 
bio-transform  under  laboratory  conditions.  We  think  that 
phototransformation  is  occurring  on  our  outdoor  samples  because  (l)  the 
IM  pieces  changed  color,  (2)  an  unknown  peak  occurs  in  the  HPLC  where 
polar  substances  elute,  and  (3)  this  peak  is  largest  in  samples  collected 
during  the  summer  months  when  sun  irradiance  is  highest.  We  expect  the 
dissolved  contaminant  mixtures  and  the  low  pH  of  water  with  NTO  to  af¬ 
fect  microbial  communities  in  ways  that  are  currently  unknown.  The  solu¬ 
bility  of  most  IM  constituents  is  higher  than  that  of  TNT  and  RDX,  in¬ 
creasing  the  likelihood  that  they  will  reach  groundwater. 

We  showed  that  both  NTO  and  DNAN  phototransform  when  in  solution. 

In  pure  water,  DNAN  phototransformed  more  than  NTO.  However,  rates 
were  sensitive  to  solution  chemistry.  NTO  phototransformation  was  en¬ 
hanced  in  acidic  and  basic  solutions.  NTO  phototransformation  rates  also 
increased  three  times  in  the  presence  of  OM,  which  is  common  in  soil  so¬ 
lutions  and  in  surface  waters.  DNAN  was  not  affected  by  OM.  We  found 
the  phototransformation  of  NTO  to  be  relatively  insensitive  to  temperature 
between  35°C  to  70°C.  On  the  other  hand,  phototransformation  rates  of 
DNAN  increased  with  increasing  temperature. 

Both  NTO  and  DNAN  experienced  adsorption  and  transformation  in  soils. 
NTO  was  weakly  adsorbed,  with  adsorption  coefficients  lower  than  those 
measured  for  RDX,  the  explosive  compound  it  is  replacing.  DNAN,  on  the 
other  hand,  was  more  strongly  adsorbed  than  TNT.  NTO  adsorption  was 
not  influenced  by  OC  in  the  soil  but  was  strongly  affected  by  soil  pH.  Ad¬ 
sorption  decreased  as  soil  pH  increased,  probably  due  to  changes  in  pro¬ 
tonation  (and  therefore  charge)  of  both  the  soil  surfaces  and  NTO.  DNAN 
adsorption  positively  correlated  with  OC  while  NTO  transformation  rates 
increased  with  increases  in  soil  OC.  Column  studies  supported  batch- 
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determined  parameters  but  indicated  that  an  increase  in  transformation 
rates  may  be  happening  in  the  soil  if  the  soil  becomes  slightly  anaerobic. 
NTO  had  little  affinity  for  soils,  making  it  more  mobile  in  the  environment 
than  DNAN;  and  as  it  is  also  more  soluble,  we  can  expect  it  to  reach 
ground  water  faster.  DNAN  is  easily  phototransformed  and  adsorbed  in 
the  soil,  making  it  less  mobile  in  the  environment. 
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